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We derive here a new stability criterion for two-fluid interfaces. This 
criterion ensures the existence of "stable" local solutions that do no break 
Qi^ ' down too fast due to Kelvin-Helmholtz instabilities. It can be seen both 

•^r , as a two-fluid generalization of the Rayleigh- Taylor criterion and as a 

nonlinear version of the Kelvin stability condition. We show that gravity 
can control the inertial effects of the shear up to frequencies that are 
C^ , high enough for the surface tension to play a relevant role. This explains 

why surface tension is a necessary condition for well-posedness while the 
(low frequency) main dynamics of interfacial waves is unaffected by it. 
In order to derive a practical version of this criterion, we work with a 
nondimensionalized version of the equations and allow for the possibility 
^ , of various asymptotic regimes, such as the shallow water limit. This limit 

being singular, we have to derive a new symbolic analysis of the Dirichlet- 
Neumann operator that includes an infinitely smoothing "tail" accounting 
for the contribution of the bottom. We then validate our criterion by 
comparison with experimental data in two important settings: air-water 
in ' interfaces and internal waves. The good agreement we observe allows us 

to discuss the scenario of wave breaking and the behavior of water-brine 
interfaces. We also show how to rigorously justify two-fluid asymptotic 
models used for applications and how to relate some of their properties to 
Kelvin-Helmholtz instabilities. 



^ ■ 1 Introduction 



1.1 General setting and overview of the results 

We are interested here in the motion of the interface between two incompressible 
fl-uids of different densities p"*" > p~ , with vorticity concentrated at the interface, 
and at rest at infinity. This problem is sometimes called the Rayleigh- Taylor 
problem; the limit cases p'^ — p~ and p~ — are also known as the Kelvin- 
Helmholtz and water waves problems respectively (see [5] for a recent review). 
The mathematical analysis of the stability issue of such interfaces has raised 
a lot of interest. The particular case of the water waves problem is certainly 
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the best understood. Following the pioneer works of [U HU IHU HH [TO], S. 
Wu established the well posedness of the water waves equation for ID and 2D 
surface waves [5S1 [57] in infinite depth, and without surface tension, provided 
that the Rayleigh- Taylor criterion is satisfied, 

(Rayleigh-Taylor) - d,P |^^^^^^^ > 0, 

where z is the vertical coordinate and P the pressure (she also established that 
this condition is satisfied as long as the wave does not self- intersect) . It is in- 
structive to remark that the linearized version of this criterion (around the rest 
state) is simply p~^g > 0, where g is the (vertical) acceleration of gravity - in 
other words, water must stands below the interface. Other approaches have 
also been developed, such as [36^ (finite depth, Eulerian framework), [39ll47l[T8] 
(fluid droplet in zero gravity, Lagrangian framework, irrotationality assumption 
removed), [TTj (weak solutions in the whole space with density discontinuity at 
the surface), etc. More recently, various works gave some additional information 
on the solution. It was shown in [31 [31] that the existence time is uniform with 
respect to the shallowness parameter (among others), allowing for the justifica- 
tion of various models used in coastal oceanography. In another direction, Wu 
[Sni [BO] and Germain-Masmoudi-Shatah [55] showed almost global {2D) and 
global (3D) existence results under smallness conditions. 

In presence of surface tension, the analysis is technically more difficult, but 
there is no striking difference for the behavior of the solution. Local well- 
posedness remains true and the solution of the water waves problem with surface 
tension converges to the solution of the water waves problem without surface 
tension when the surface tension coefficient goes to zero [62l [S] [6] [40]. It is 
not known however whether this zero surface tension limit is compatible with 
the shallow water limit; the justification proposed in [351 [301 [31] for instance of 
some shallow water models in presence of surface tension is relevant for large 
values of the surface tension (more precisely, for Bond numbers of order 0(1)) 
but provide an existence time extremely short when the Bond number is large 
(small surface tension) , which is the case in coastal oceanography for instance. 

For the Rayleigh-Taylor problem (p^ > 0) the picture is completely different. 
It is known, at least for ID interfaces, that, outside the analytic framework 
of [55] [54] , the evolution equations are ill-posed in absence of surface tension 
[23} l33l l34] - see also [38l [58] for the related Kelvin- Helmholtz problem where 
yO"*" = p~ - even in very low regularity (C""^"*"" in |34| I38j and chord-arc in 
[58]). The reason of this ill-posedness is that the nonlinearity creates locally 
a discontinuity of the tangential velocity at the interface that induces Kelvin- 
Helmholtz instabilities. 

More recently, it has been proved that taking into account the surface tension 
restores the local well-posedness of the equations [H [7] [3H1 US]- However, the 
existence time of the solution provided by these results is very small when the 
surface tension is small. The fact that the role of gravity (or gravity itself) 
is not considered in these references suggests that these general results can be 
improved in the "stable" configuration where the heavier fiuid is placed below 



the lighter one. In particular, we would like to explain two physical phenomena 

t|l Air-water interface. In coastal areas, for instance, waves propagate over 
several wavelength (depending on various physical parameters introduced 
below). This observation is consistent with the results obtained in [3] 
where the air density is neglected {p~ =0). If this density, which is small 
but non zero, is taken into account, the existence results of [4j [71 l48l fT3] 
ensure that surface tension can prevent Kelvin-Helmholtz instabilities, but 
for times that allow the wave to propagate over a few millimeters only. 
We therefore want to understand the mechanism that controls Kelvin- 
Helmholtz instabilities for time scales consistent with the observations. 

t|2 Internal waves. Waves at the interface of two immiscible fluids of pos- 
sibly very similar densities have been widely investigated in the physics 
literature. Waves at the interface of two layers of water of different den- 
sities in the ocean are also commonly described with this formalism. As 
for the case of the air-water interface above, the known results only pro- 
vide an existence time smaller by many orders of magnitude than what is 
observed. 

These two phenomena raise an apparent paradox: 

P The asymptotic models describing waves in coastal areas or internal waves, 
and which provide very accurate results, all neglect surface tension. On 
the other hand, the analysis of the full equations tells us that surface 
tension is crucial for the existence of a solution. In other words, surface 
tension allows interfacial waves to exist but does not play any role in their 
dynamics. 

Very recently, some progress has been made to understand jjl. It has been proved 
in [m |45i that, the surface tension coefficient being kept fixed, the solution of 
the two-fluid system converges to the solution of the water waves problem (with 
surface tension) as the density p~ goes to zero. In [iS], the following important 
result was proved: the solution of the two-fluid system converges to the solution 
of the water waves equations (without surface tension) when both the density p~ 
and the surface tension go to zero. It is just required that (p^)^ < a'^^'^, where 
a is the surface tension coefficient. A corollary of the result presented here is 
a slight improvement of this condition, showing that the result holds under the 
weaker condition (p~)^ = o(cr) as suggested by the linear analysis. However, in 
physical applications, the ratio (p~)^/(T cannot be made arbitrarily small, and 
it would be interesting to have at our disposal a criterion to check whether a 
given interfacial wave will be "stable" (in the sense of "observable on a relevant 
time scale" ) or will almost instantaneously break down due to Kelvin-Helmholtz 
instabilities. 

The need for such a criterion is also crucial to explain tJ2. Indeed, the small- 
ness of the ratio (p^)'^ /a obviously cannot be invoked to explain the observation 
of such phenomena as internal waves since the relative density of the upper fluid 
is not small in that case (while a remains very small) , and a better understand- 
ing of the formation of Kelvin-Helmholtz instabilities is required. 



A hint of what such a criterion should look like can be obtained by consid- 
ering the related problem of the linear stability of two fluid layers moving with 
constant horizontal velocities c"*" 7^ c~ (with c^ € M'', d = 1, 2). It is known 
since Kelvin that small perturbations of the interface are linearly stable if the 
following criterion is satisfied (see for instance [T^ ) 

(Kelvin) g{p+ p-) > ^ Jf^^c(0)|Ic±]|4, 

where [c ] = c^ — c'' and c(0) is a constant taking into account the geometry 
of the problenu- The equivalent of this linear criterion in the water waves case 
is the linear Rayleigh- Taylor criterion already mentioned, namely, p'^g > 0. In 
the problem we investigate here both fiuids are at rest at infinity, but there is is 
a local discontinuity of the horizontal velocity created by the nonlinear motion 
of the interfacial waves. We can therefore expect that a local version of the 
Kelvin criterion may be derived to assess the local linear stability of interface 
perturbations. A natural question is then: is there a nonlinear generalization 
of this expected local Kelvin criterion ? Or equivalently, is there a generalization 
to two- fluids interfaces of the Rayleigh- Taylor criterion? 

The first main result of this paper is to show that such a nonlinear criterion 
exists. It can be stated as 



1 {p+p-f 



9^P^ U^cl > T^V-^TTo c(C) HL, (1-1) 



where ^ is the interface parametrization, u — fV^ \z=A ^^ ^^^ jump of the 
horizontal velocity at the interface, and c(C) is a constant that depends on the 
geometry of the problem (two layers of finite depth in this paper) and that can 
be estimated quite preciseljo. 

The question is now to check whether the criterion (jl.ip can explain (for 
instance) the observation of the phenomena ttl and t|2 mentioned above. 
For ttl (air- water interface), it is likely that the relative density p~ = p~ /[p'^ -\- 
p~) is small enough to make the r.h.s. of (jl.ip smaller than the l.h.s. even 
though a is small. This is at least the case if we consider the limit {p~Y /a — > 
as in [i^ . 

For tt2 (internal waves), the relative density p~ is close to 1/2 and cannot be used 
to explain why (|1.1|) should be satisfied in that case. Therefore, the smallness of 
the discontinuity on the horizontal velocity is the only reason that could explain 
tt2. 

The main idea behind the stability criterion (jl.ip is the following expla- 
nation of the apparent paradox p, that we can state in rough mathematical 
terms as follows: the Kelvin- Helmholtz instabilities appear above the frequency 
threshold for which surface tension is relevant, while the main (observable) part 



^c{0) = 1 for the case of two layers of infinite deptfi considered in | 12 |: see Remark 1 5. II for 
an explicit expression in the case of two layers of finite depth 

^for a flat interface C = 0, this constant is the same as in the linear Kelvin criterion 



of the wave involves low frequencies located below this frequency threshold; con- 
sequently, the Kelvin-Helmholtz instabilities are regularized by surface tension, 
while the main part of the wave is unaffected by it. The role of criterion (jl.ip is 
to ensure that the frequency threshold mentioned above is high enough. For low 
frequencies, we show that gravity can stabilize the inertial effects of the shear. 
The main task in this paper is to prove rigorously this scenario. 

Remark 1.1 (Rayleigh- Taylor stability vs. Kelvin-Helmholtz instability) . As al- 
ready mentioned, surface tension has been shown to regularize Kelvin-Helmholtz 
instabilities (that are due to the shear). As a matter of fact, it also stabilizes 
Rayleigh- Taylor instability in the sense that it allows for a local solution even if 
the heavier fluid is placed above the lighter one. Such a solution is an example 
of what would be an unstable solution for us: even though the initial value 
problem is locally well-posed, it does not satisfy (|1.1|) and the existence time is 
very small. In the stable configuration (heavier fluid below), we said that we use 
gravity to control the destabilizing effects of the shear; equivalently, Rayleigh- 
Taylor stability controls the Kelvin-Helmholtz instabilities at low frequencies. 

Having proved that (jl.l[) determines the stability of interfacial waves (and 
constructed such solutions to the two- fluid equations), one practical question 
arises: how can we have access to the quantities |— 9zP^] and |w|oo that we 
need to know to determine whether (ll.ip is satisfied or not? Indeed, one would 
like to use this criterion to assess the stability of some given configuration: for 
instance, knowing the depth and density of both fluid layers, is it possible to 
have a stable perturbation of the interface of amplitude a and wavelength A? 
It is quite difhcult to know |— cJ^P^I and |a;|oo experimentally (see however [27] 
for measurements oiuj). A useful approach is therefore to perform an asymptotic 
analysis of the problem. Of particular importance for applications is the shallow 
water regime (the wave length of the perturbation is large compared to the 
depths H~^ and H^ of the fluid layers) . In this regime, the pressure satisfies the 
hydrostatic approximation with a fairly good precision, and one has therefore 
|— c^^P*] ^ (p"*" — p~)g, it is also possible to show that w has a typical size 
of order jj-^g'H (with g' = (/?+ — p^)g and H — +H-+n-H+ )' ^ prediction 
that is consistent with the measurements of [37] for instance. Plugging these 
approximations into p.ip shows the relevance of the dimensionless parameter 
T defined as 



_^^a^ {p++p~)g' 



and suggests a very simple practical stability criterion in shallow water 

T <C 1 : Stable configuration; T ^ 1 : Unstable configuration. (1-2) 

When T ~ 1, stability/instability is critical and it is necessary to look at the 
exact criterion (11.11) . 



In order to prove rigorously the relevance of the practical criterion (|1.2p . 
one must be able to handle the shallow water limit in the construction of the 
stable solutions allowed by (II. ip . Unfortunately, this limit is singular, which 



complicates the proof. In particular, it restricts the number of tools at our dis- 
posal since standard symbolic analysis is not adapted to shallow water regimes. 
Roughly speaking, symbolic analysis neglects the information coming from the 
bottom since it is infinitely smoothing (this is the argument used in [T] to prove 
local well posedness of the water waves equations over very exotic bottoms). 
However, in the shallow limit, the influence of the bottom is the main factor in 
the evolution of the wave: though infinitely smooth, its contribution becomes 
quantitatively very large. In the water waves case, various techniques allow one 
to bypass symbolic analysis and justify some shallow water models [21 [31] . In the 
two-fluid case under investigation here, it is not clear whether these techniques 
can be adapted. We therefore appealed to a technical result of independent 
interest, namely, a symbolic analysis "with tail" of the Dirichlet-Neumann op- 
erator, the tail corresponding to the infinitely smoothing component of this 
operator taking into account the influence of the bottom. This allows us to use 
some (rudimentary) symbolic analysis without suffering from the shallow water 
singularity. 

Once this analysis is done, we know that interfacial wave are well behaved 
in the shallow water limit when the stability criterion (jl.ll) — or its practi- 
cal version (|1.2p — is satisfied. We then turn to apply this result to the two 
phenomena (tl and (12 described above. The good agreement with experimen- 
tal data allow us to confirm the relevance of one-fluid asymptotic models for 
air-water interfaces, and to give the first rigorous justification (on the relevant 
time scale) of two-fluid asymptotic models commonly used to describe internal 
waves. We also indicate how some singularities of these models can be related 
to Kelvin-Helmholtz instabilities. We finally discuss two phenomena: 

t|l Air- Water interfaces. Simple physical examples show that close to wave 
breaking, it is likely that p.ip fails to be satisfied for certain configura- 
tions. This would mean that the first singularity observed in wave breaking 
may not always be a singularity of the water waves equations, but some- 
times a two-fluid singularity. 

tt2 The description of water-brine interfaces (as for oceanic internal waves) 
does not fall into the range of equations (|1.3p - (|1.8l) because water and 
brine are not immiscible. We use the stability criterion(s) derived in this 
paper to propose a two-fluid description of the water-brine interface that 
includes an artificial surface tension a. Using comparison with laboratory 
experiments we are able to propose a numerical value for a. 

In order to make a more detailed description of the results described above, 
we first need to write the equations, and to derive their nondimensionalized 
version. This is done in the following two subsections. 

1.2 The equations 

We assume that the interface is parametrized by a function ({t,X) {X e R'') 
and denote by fij and fi^ the volume occupied by the lower and upper fluids 



respectively at time t. Choosing the origin of the vertical axis to correspond 
with the interface between the two fluids at rest, we assume that ilf (resp. $7^) 
is bounded below (resp. above) by an horizontal wall located at z — — iJ+ (resp. 
z = H~). We also denote by Tt the interface Vt = {{X, z), z — ({t, X)} and by 
r the upper and lower boundaries T = {z = ^H^}. 
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Finally, we denote by U the velocity field in J7( ; the horizontal component 
of U* is written V^ and its vertical one w^ . The pressure is denoted by P^ . 
For the sake of clarity, it is also convenient to introduce some notation to express 
the difference and average of these quantity across the interface. 

Notation 1.2. If A~^ and A^ are two quantities (real numbers, functions, etc.), 
the notations |A*] and (A^) stand for 

lA^j = A+-A- and (A±) = ^^ + ^^ . 

We can now state the equations of motion. 

• Equations in the fluid layers. In both fluid layers, the velocity field U* 
and the pressure P* satisfy the equations 

div U±(i,-) ==0, curl U±(t, •) =0, infl^ (* > 0), (1.3) 

which express the incompressibility and irrotationality assumptions, and 

p±(atU± + (U±.Vx,.)U±)=-Vx,.P±-.9e, inftf (t > 0), (1.4) 

which expresses the conservation of momentum (Euler equation) . 

• Boundary conditions at the rigid bottom and lid. Impermeability of these 
two boundaries is classically rendered by 



«;±(t,-)|,± =0, 



{t > 0). 



(1.5) 



• Boundary conditions at the moving interface. The fact that the interface is 
a bounding surface (the fluid particles do not cross it) yields the equations 

dtC - Vl + |VCPU± - 0, (i>0), (1.6) 

where U^ :— U^ij, • n, n being the upward unit normal vector to the 
interface Ft. A direct consequence of (II. 6p is that there is no jump of the 
normal component of the velocity at the interface. Finally, the continuity 
of the stress tensor at the interface gives in our particular case 

[P±(i,.)|,J-c7^(C), (i>0), (1.7) 

where a is the surface tension coefficient and ^(C) denotes the mean cur- 
vature of the interface, 



It is quite convenient for the mathematical analysis of the equations (jl.3p - 
(|1.7I) to transform them into a set of scalar equations on the interface. From the 
irrotationality assumption, we can write the velocities U in terms of velocity 
potentials $*, namely, 

U±(i,-)=Vx,.$*(i,-) i^^f (i>0). 

Defining the trace of $* at the interface by 

V'±(i, •) = <&*(*, Olr, (i>0), 

we can now reduce the equations (jl.3p - (|1.7p to a set of equations on C, and 4' , 
namely, 

dtC-g+m+^0, (1.8) 
g + [(:]ij+^g-m-, (1.9) 

9C + V^^P-^^^i^i^^^^^^,?^)^-P^ (1.10) 

■^^2' ^ ' 2(1 + |VCP) / kt' ^ ' 

[P±(f,.)|rJ = ^^(C), (1.11) 



P^{dt^^ 



where § * [C] are the Dirichlet- Neumann operators corresponding to the two 
fluid layers and defined (under reasonable assumptions on C and i/j^) by 



c^±[C]^± = Vi + |vcpa„$±|,^, 

where $* solve 



Ax,.$* = in nf, 






and 9„$| stands for the upward partial derivative of $* at the interface. 



We now show|j that it is possible to reduce the two-fluid equations to a set 
of two equations on the surface elevation (^ and of the quantity ip defined as 

"0 := p^ijj^ — p^ip~ , 

where p^ stands for the relative density, p^ = 4,^ _ (in particular, p^+p^ = 1 
and p^ — p^ is the Atwood number) . Defining the operator ^ [C] by 

one can write (^ +[(']?/'+ — g^[(]ip^ in terms of C and ip only, 

and one can also get ip^ in terms of C and ip, 

i^^ = g^[Cr'og[C]i,. (1.15) 

A formulation of the two-fluid equations as a system of two scalar evolution 
equations on C and ip can then be given. 



r dtc-g m = 0, 



dt^ + g'C + 2^*1^^* 



i[p±((?±[C]V±+vc-v^±)2] 



-^(C), 



2 1 + lVCP P++p- 

(1.16) 
where ip"^ = g'^lQ ^ o g[(]ip and g' stands for the reduced gravity, 

9' = {P" - pr)9- 

Remark 1.3. In the case p~ = 0, one has "0 = V'^ and g [Qip = g^[Qi^^ and 
(|1.16l) reduces to the formulation of the water waves equations in terms of the 
surface elevation C, and the velocity potential at the surface ip^ due to Zakharov 
[55] and Craig- Sulem [21]. In the case p~ ^ Q, C, and '0 are the two canonical 
variables of the Hamiltonian formulation exhibited by Benjamin and Bridges [5] 
(see also [ID]) 

1.3 Nondimensionalized equations 

The order of magnitude of the existence time of the solution is our main concern 
in this paper. It is therefore crucial to get some information on the size of the 
different quantities appearing in the equation, and in particular of those that 
play a role in our stability criterion (|l.ip . These informations can be obtained 
by experimental measurements, but it is of course preferable to now them a 
priori in terms of the physical characteristics of the flow under consideration. 



^We do not justify the different steps of the derivation here; this will be done in Section [2] 



namely, 

- The typical size a of the interfacial waves 

- The typical wavelength A of these waves 

- The depth H^ of the upper and lower fluid layers at rest 

- The value of the surface tension coefficient a 

- The densities p^ of the two fluids 

We show in Appendix \X\ how the linear theory can be used to obtain such an 
a priori estimate of the order of magnitude of the different unknowns. The best 
way to exploit this information is to nondimensionalize the equations, that is, to 
perform a linear change of unknowns and variables such that all the unknowns 
are now dimensionless and of size 0{1) for typical configurations. 

Writing (J1.16p in such a dimensionless form requires the introduction of 
various parameters; the first two are given by 

a iJ2 H+H- 

with H 



H' ^ A2 ' P+H-+P-H+' 

in the water waves configuration, i.e. when p~ = 0, e is called the amplitude or 
nonlinearity parameter and /i the shallowness parameter. Throughout this arti- 
cle, fi is assumed to remain bounded, as well as s. For notational convenience, 
and without loss of generality, we take 

0<Ai<l, 0<£<1. 

The above definition of e and /i takes into account both fluid layers at the same 
time. For a more specific description of each fiuid layer, it is convenient to 
introduce 

^ " i/± ' ^ A2 • 

We also need to define the relative depth H and the Bond numbeQ Bo as 

— H ' a 

Remark 1.4. In full generality, one should allow H'^ to be much larger or much 
smaller than H~ (as in [11] for instance). For the sake of clarity, we assume 
throughout this paper that _ff + and H~ are of same order (so that e+ ~ e~ ~ e 
and fi'^ ~ fi~ ~ /i). There does not seem to be any obstruction other than 
technical to generalize our method to more general configurations. 

The next step is to introduce dimensionless Dirichlet-Neumann operators. 



g^leCK^]^^ = v/i + |evcpa„$±^^^. 



*The Bond number measures the ratio of gravity forces over capillary forces. For some 
reason, it often appears as tlie inverse of this quantity in the mathematics literature 
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where $^ solve 

and 9n^^ i^^^;; stands for the upward conormal derivative of <l>^ at the interface 
(that is, g + [eQil)^ = dz^^ |^^^^ -e^VC- V(/)=^ \.=ec:)- The operator Q [C] defined 
in (|1.13l) has similarly a dimensionless version given by 



gM]^<^M,H-]{p+g-[eC.H-]-p-g + [eQ,H+]] g + [eCH+]. (1.18) 



We show in Appendix |^ that the equations (J1.16p can then be written in 
dimensionless form as 

(1.19) 

with V'± = g+H,H+]-^ o ^^KIV'- 

1.4 Description of the results and organization of the pa- 
per 

We give in Section [2] some preliminary results that give a rigorous basis to 
the manipulations performed in ? ^1.2I and 11.31 above to derive the equations 
(|1.16p and (J1.19p . General notations and definitions are given in J|2.1l and the 
Dirichlet-Neumann operators g^ [eC,, H ] are studied in Ji2.2l Their invertibility 
properties are investigated in f H2.2.3p : a technical difficulty is that the inverse 
of the DN operator takes values in Beppo-Levi spaces (i.e. functions with gra- 
dient in Sobolev spaces, but not necessarily in L^ since we consider unbounded 
domains). The structure of the shape derivatives of the DN operators are then 
discussed in i )2.2.4l We are then able to address in i j2.3l a transmission problem 
that is crucial to derive rigorously the two-fiuid equations p.l9p since it states 
that the trace of the potential at the upper and lower part of the interface, i/'" 
and ip'^ , are fully determined by C and "ip = p'^ip^ — p^ip"- This allows us 
to give a rigorous definition to the operator g — gfj.[eC] introduced in (jl.lSp . 
Since we want to be able to handle the shallow water limit, we need to give two 
versions for (almost) all the estimates given in this section: one that is sharp 
with respect to regularity but not optimal with respect to the dependence on 
/i, and one that is sharp with respect to p but not optimal with respect to 
regularity. For instance, gti^CTlL ] has an operator norm of size 0{y^) if it is 
seen, as usual, as a first order operator; but if we consider it as a second order 
operator, we get a better control of size 0{fj,) on its operator norm (see Remark 
O below). 
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Section [3] is devoted to the symbolic analysis of the Dirichlet-Neumann op- 
erator g+ = H'^g;i:[eC,H+]. it is standard [H [36l |47l [2] that the principal 
symbol oi Q^ can be written in terms of the Laplace-Beltrami associated to the 
surface. More precisely, if we define the symbol g{X,£,) as 



then the following holds, 

VO<s<io, \g+i'~^gix,D)^\H.+i/2<M{to + 3)\\7^\H.-i/2, (1.20) 

where to > d/2 and M{to + 3) is a constant depending on |C|_f/*o+3 and the 
minimal depth of the lower fluid. Since Q "*" and g{x, D) are first order operator, 
this identity shows that Q^ can be replaced by g{x,D) up to a more regular 
(zero order) operator. Unfortunately, such a substitution is not possible in the 
shallow water limit since this zero order operator is of size 0(1) with respect 
to /i, while both g^ and g{x,D) have operator norms of size 0{y/JI). Using 
symbolic analysis induces therefore a singularity of order 0(/x~^/^) when the 
shallow water limit is considered. The explanation of this behavior is that the 
shallow water regime corresponds to configurations where the bottom plays a 
very important role; however, the contribution of the bottom to the Dirichlet- 
Neumann operator is analytic and therefore neglected (at any order) by symbolic 
analysis. In order to handle this difficulty, we propose a symbolic analysis "with 
tail" of the DN operator, that takes into account the infinitely smoothing (but 
very large) contribution of the bottom. This leads us to replace the symbol 
ff(X,0 by ff(X,Ot+(X,e), where 



t+(x n - a + .+n /" VieP + ^M^ + i)^(lvcn^P-(vc-e)V 

t {X,0-{l + e 0J_^ l + £2^(z + l)2|VCP ^' 

(when d = 1, this simplifies into t+(x,n = (1 + e+O ^'''^ '^"^ \^\)- We 
then show in Theorem 13.11 that (|1.20p can be improved into 

'^^'""'"' -1/2. 



I^+V- VA^Op[5tanh(VM+i+)]V'|jj.+i/2 < £VJ^M{to + S)\ViP\h. 

The cost in terms of derivatives of tp is the same as in (|1.20p , but the behavior 
with respect to the parameters e and fi is much better since there is a gain 
of size 0{ey/JI). This allows us to handle the shallow water limit. We also 
give in Section |3] the symbolic analysis (with tail) oi {g~)~^g'^, {g~)~^g and 
Q {g^)^^dj that are needed for the analysis of the two-fiuid equations. 

In Section m we show that it is possible to quasilinearize the two-fiuid equa- 
tions (I1.19|) by differentiating them and writing them in terms of the relevant 
unknowns. We first introduce and study in i j4.1l some new operators that are 
needed to write the quasilinearized equations. An important step is a lineariza- 
tion formula for the operator Q that is given in ij4.2l the main step to establish 
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this linearization formula is an explicit formula for the derivatives of the map- 
ping ^ i-7> gi^[eQtjj (shape derivative). This linearization formula also suggests 
what the good unknown should be. More precisely, the equations (J1.19p differ- 
entiated a-times should not be written in terms of 9"C and d°'^j but in terms 
of C(q) and -0(0) defined as 

where w is the vertical component of the velocity at the interface in the ± 
fluid, while ip^ is as in ()1.19p . When p~ = 0, these unknowns coincide of course 
with the ones used in 32, 52i to write the quasilinearized water waves equations. 
The quasilinearized equations are then derived in §4.31 Differentiating a times 
the second equation of (|1.19p and rewriting it in terms of C(q) and V'(a) , one gets 
(without surface tension) 

at0(„) + ad^C + siV"^) ■ VV(a) + el^^l • V(p±Vf„)) - 

(the symbol ^ means that harmless terms are omitted); in this expression, a 
can be related to the (dimensionless) jump of the vertical derivative of the pres- 
sure at the interface, and V_ stands for the horizontal velocity at the interface 
in the fluid ±. In order to write this equation in terms of the good unknowns 
C(q) and V'(q)j it is necessary to write {p^ipfa)) ^^ terms of Ip^V'fo.)! ~ ''/'(") ^-'^'^ 
C,{a)- This is completely trivial in the water waves case p^ = 0. The depen- 
dence of this term on C,ia.) is therefore specific to the two-fluid system, and it is 
responsible for the Kelvin- Helmholtz instabilities — through the operator £ [C] 
in the quasilinear system given in Proposition 14.91 

The main existence and stability results are then given in Section [S] We first 
state in i j5.ll several versions of the stability criterion. The first one (see ijS.l.ip 
differs slightly from (|l.ip because it also involves the L°°-norm of first order 
space-time derivatives of the jump of the horizontal velocity. In §5.1.21 we show 
how to recover (jl.ip when the jump of velocity is nonzero. The practical criterion 
(|1.2p is then derived in ^15.1.31 After some considerations on the conditions 
wc impose on the initial data (see H5.2p . we state in H5.3I the main results of 
this paper, which show that "stable" solutions exist under the criterion given 
in § §5. J.JI and \5.lM (we give two versions of the theorem). By stable, we 
mean that the existence time depends on the surface tension a through the 
stability criterion (|5.1|) only. In particular, it is not necessarily small if a is 
small. We want to emphasize on the fact that the theorems are given under a 
nondimensionalized form and that the existence time is uniform with respect 
to the parameters e and /i. In particular, this allows one to handle the shallow 
water limit (or other regimes) straightforwardly. Other kinds of limits (zero 
density and zero surface tension) are also direct corollaries. We then prove in 
^5.41 that these stable solutions persist over large times (i.e. of size 0(l/e)) if a 
stronger version of the stability criterions is satisfied. 

The proofs of the theorems are then given in ^15. 51 The stability criterion appears 
in the key Lemma 15.111 where it is used to control the destabilizing effects of 
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the shear by the surface tension and gravity terms. As in the water waves case, 
some technical difficulties arise because of the second order surface tension term 
that requires special care in the treatment of the subprincipal terms in the 
symmetrization process. We chose to adapt the technique introduced in [52] 
(see also [S3] for a synthetic presentation), which is probably not the sharpest 
one in terms of regularity (see for instance [JT] [S] [HI UHl HI HI] for alternatives 
methods) but seems to be the most robust with respect to the shallow water 
limit, in particular because it does not require symbolic analysis. This technique 
requires that time derivative be treated as space derivatives, and this is the 
reason why this is done throughout this paper. 

We finally sketch in Section [5] some applications for the stability criterion 
(jl.ip and its practical version (11.21) . We consider in Ji6.1l the case of air- water 
interfaces characterized by a very small density of the upper fluid. We first 
confirm in H6.1.1I that the density of the air can be neglected in the asymp- 
totic models used in coastal oceanography to describe the propagation of waves 
(we consider here the example of a so called long wave). For waves close to 
the breaking point, we show in i)6. 1.21 that this might not be true anymore and 
that Kelvin-Helmholtz instabilities may appear. We therefore suggest that wave 
breaking may sometimes be a two-fluid singularity rather than a singularity of 
the water waves (one-fluid) equations. Links with physical phenomena such as 
spilling breakers and white caps are also commented. 

We then consider in t ^6.2l the case of internal waves at the interface of two fluids 
of comparable density (stability of the fluid is therefore ensured by a small shear 
velocity). We first check that some "stable" configurations reported in exper- 
imantal works do satisfy our stability criterion. We then consider the case of 
water-brine interfaces, of interest for the study of oceanographic internal waves. 
There is no natural value for the surface tension for such interface (water and 
brine are not immiscible) , but investigating the occurrence of Kelvin-Helmholtz 
instabilities thanks to the experiments of [2 7) . we show how to use our stabiliy 
criterion to propose a value for the artificial surface tension that must be intro- 
duced to use the two-fluid formalism. 

We end this section by showing how to justify rigorously the two-fluid asymp- 
totic models used in applications, and on a time scale consitent with the observa- 
tions. We give the details for a shallow-water/shallow-water model and explain 
how one of its singularities can be related to Kelvin-Helmholtz instabilities. 

1.5 Basic notations 

- We denote by C(Ai , A2 , . . . ) a constant depending on the parameters Ai , A2, . . . 
and whose dependence on the Xj is always assumed to be nondecreasing. 

- We denote by X = {Xi, . . . , Xd) € R'' the horizontal variables, by z the verti- 
cal one, and by t the time variable. 

- We denote by dj {1 < j < d) partial differentiation with respect to Xj; partial 
differentiation with respect to t and z is denoted by dt and dz. 

- We denote by V and A the standard gradient and Laplace operators in the 
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horizontal coordinates. These operators are denoted by V x.z and Ax,2 when 
they also take into account the vertical variable. 

- We denote by V^ = {y/JIV'^ , dz)"^ , that is V^ corresponds to Vx.z with a 
factor y//i in front of each horizontal derivative. 

- We denote A := (1 - A)!/^ and H''{R'^) {s e R) the usual Sobolev space 
H^R'') = {u e S'{R'^),\u\h^ < oo}, where \u\hs = \A'u\l2. We keep this 
notation if u is a vector or matrix with coefficients in H''{M.'^). 

- We denote by (•, •) the standard L^-scalar product. 

- The notation a V 6 stands for max{a, b}. 

- We denote by Op(cr) or a{x,D the pseudodifferential operator of symbol 



2 Preliminary results 

2.1 Notations and definitions 

2.1.1 General notations 

Throughout this section, we always assume that the interface deformation ( G 
H*°^^(W^) (to > d/2) does not touch the bottom nor the lid. In dimensionless 
variables, this condition reads 



3h^^„ > 0, inf (1 ± e^aX)) > h^^^. (2.1) 



We also denote hy fl the domains 

n+ = {{X,z)eR'^+\-l<z<e+C{X)} 
n- = {{X, z) e M''+\ e-C(^) < 2 < 1}. 

We finally introduce, for notational convenience, a constant M defined aCl 



1 iJ+ H- 



M-C(^^,— ,— ,|CU.o+2), (2.2) 



as well as 

M(s) = C(M,|C|hO, (2.3) 

where we recall that C(-) denotes generically a nondecreasing, positive, function 
of its arguments. 



^the dependence on H^ /H~ and H~ /H~^ is harmless because we assumed that _ff+ and 
H~ are of same order (see Remark 1 1.411 . For more general configurations, a finer analysis on 
the dependence on H^ of the solution is needed 
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2.1.2 DifFeomorphisms 

It is often convenient to transform boundary value problems on the fluid domains 
n^ into boundary value problems on the flat strips S^ defined as 

S+ = {iX,z),-l< z<0} and 5~ = {(X, z),0 < z < 1}, 

so that S^ correspond do the two fluid domains at rest. Various diffeomorphisms 
can be used for such domain transformations. We only use here admissible 
diffeomorphisms in the following sense: 

Definition 2.1 (Admissible diffeomorphisms). Let to > d/2 and C e i7*»+2(R'^) 

be such that (|2.II) is satisfied. 

We say that E* : S^ — > 57* is an admissible diffeom,orphism, if 

1. S* can be extended to the boundaries in such a way that 

E±({z-0}) = {z = e±C} and S±({z = t1}) = {C = Tl}- 

2. The coefficients of the Jacobian matrix J^i are bounded on S^ , and 

I^S±Il~(5±) <M. 

3. The determinant | det J^i | is uniformly bounded from below on S^ by a 
nonnegative constant c* such that 4r < M. 

If $* is defined on ft^, then one can define </)* — $* o S* on S^; a direct 
application of the chain rule shows that 

where /^ is the (d + I) x {d + I) diagonal matrix with entries \/fJ^ on the 
d-first diagonal coefficients, and 1 on the last one. 

In particular, the (variational formulation of the) equation (d^ + /i^ A)$* = 
in 51^ is transformed by E* into the (variational formulation of the) following 
variable coefficient elliptic equation on S^ , 

V^* .P(E±)V^*0=^ =0 in 5=^, (2.4) 

where 

P(S±) ^\detJ^^\iI^^)-'j-l{I'^^r{J-lf{I^^)-'; (2.5) 

it is easy to check (Prop. 2.3 of [3]) that the matrix P(E^) is uniformly coercive 
with coercivity constant fc(E^). We assume further that 

||P(S±)||oo < Af and ^^<A/. (2.6) 

It is also convenient to assume that P(E*) is a Sobolev perturbation of the 
identity matrix (this is in general a direct consequence of the fact that C has 
Sobolev regularity) in the sense that 

||A*°+i(Id -P(E±))L^^.^<M. (2.7) 
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Let us also mention that we use the notation 9„ for the transformed equations 
to denote the upward conormal derivative naturaUy associated to any elhptic 
equation. For instance, the upward conormal derivative associated to (|2.4p is 

9„</>± =e, •P(I]±)V''*0±. (2.8) 

We now introduce as in |36l (see also TO') a very useful class of admissible 
diffeomorphisms called regularizing diffeomorphisms. 

Definition 2.2 (Regularizing diffeomorphisms). Let ip > d/'^, s > and C, S 
jjt„+2 p //^s(Rd) ^^ g^pj^ ^j^j^^ (gjj -g satisfied. 

We say that an admissible diffeomorphism S^ : S^ — > $7^ is regularizing if it is 
of the form S±(X, z) = (X, z + cr±(X, z)) with 

||A^V^a±|U2(5±) + ||A^-i9,W±|U2(5±) < A/(s + 1/2). (2.9) 

Note that if S* is regularizing, then one automatically has from (j2.5p that 

||A^(F(I]±) - Id )\\ms±) < M{s + 1/2). (2.10) 

Example 2.3. The most simple example of admissible diffeomorphism is given 

by 

E±(X,z) = (X,e±(l±z)C(X)+z). 

It satisfies (EH) and dUT]), but not ^^ (the r.h.s. should be M{s + 1) in this 
inequality) . 

Example 2.4. For small enough S > (and with x : K — > R a smooth, positive, 
compactly supported function such that x(0) = 1)? the functions S defined as 

Y.^{X,z)^{X,z + a^{X,z)) with a^{-,z) ^ e^{l ± z)x{6z\D\)C 

are regularizing diffeomorphisms (see [36]) also satisfying (|2.6p and (|2.7p . 

2.1.3 Functional spaces 

We introduce here two functional spaces that are closely related to the energy 
norm. The space H^^^(R'^) is introduced to measure the extra-control (with 
respect to the iJ^-norm) provided by the surface tension term and is defined 
for iV G N as 

H^+\R'^) ^{ue H^{R'^), -^Vm e H^{Ry}, (2.11) 

vBo 

endowed with the norm |w|?^n+i = |u|^n + ";^|Vw|^n. 

Since the velocity potential does not necessary decay at infinity (while the veloc- 
ity does), we have to work in Beppo-Levi spaces [55] on open subsets il C R^^^^ 
defined for iV > as 

H^+\n) = {ue Ll^m, Vx,zu e F^(f])'^+i}, 
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endowed with the (semi)-norm |1u|1^n+i = ||Vu||/fiv. We also define similar 
spaces over R'^, 

jjs+i/2^^d^ = {m e Ll^{M.'^).yu e H'-^/^{Ry} {s e M), 

endowed with the (semi)-nonn |u|^s+i/2 = |VM|/fs-i/2. In order to capture the 
shallow-water dynamics, it is convenient to introduce a variant of this space 
that depends on the shallowness parameter fj,, 

m+^/'^{R'')^H'+'^/^{R'') endowed with |mL.+i/2 = I^JuIj^., (2.12) 
where *p is the nonhoniogeneous Fourier multipher of order 1/2 defined as 

2.2 A few results on Dirichlet-Neumann operators 

2.2.1 Two equivalent definitions 

If C e iJ*n+2(R'^) {to > d/2) satisfies ^^ then it is well known that for all 
ij}^ g iJ^/^(]R''), there exists a unique solution $='= e H'^{rt^) to the elliptic 
boundary value problem 

/ (af + /.±A)<I>±=0 in f]±, 

|*L,^=^^ d.<^U-..=0 (2.14) 

(since the interface is a graph, it is possible to view the trace of $* on the 
interface as a function defined on R'') so that the following definition of the 
Dirichlet-Neumann operators ^\ [e^Cj 1] makes sense. 






^uil*^ C4] : ,/,± ^ . /T~n7±T77l2;5 (D± 



Remark 2.5. The operators gf [eC, i? ] that appear in the introduction are 



where $^ solves p.l4p and 9„ denotes the upwards conormal derivative 

related to ^ ± [^^Ci 1] through the scaling law 

1 



g^[eCH^] = ^g^Ae^Cl]. (2.15) 



Notation 2.6. For the sake of simplicity, we simply write 



g^^g^,[s^Ci] 



when no confusion is possible. 
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From the discussion in N2.1.21 Q tp can be equivalently defined as 

^±7^±=9„(/.±l_„ (2.16) 

where (9„ now stands for the upwards conormal derivative associated to the 
elhptic operator in the fiat strip (see (|2.8p ) solved by 0^, 



V^* •F(S±)V^*0± =0 in 5± 



±. -,/,± « ^±. -n (2-1'^) 



(E* being an admissible diffeomorphism in the sense of Definition I2.ip . 

2.2.2 Basic estimates 

Let to > d/'^ and C G H*-°^'^{W^). If 1]='= is an admissible diffeomorphism in 
the sense of Definition 12.11 then for all -0^ G H'^/'^iW^), there exists a unique 
solution 0± e H'^{S^) to ([2Tf| . Moreover, one has (see [36 1 132 ] ). 

VO < s < io + 1, ||A^V^*0±|U2(5±) < AfV^|V±|^.+i/2, (2.18) 

where M is as in (|2.2p and ij^ ' (R'^) as in (|2.12p . The reverse inequality also 
holds (see Appendix IB. l.ip . 

V7I|^±|^.+i/. < M|1A^V^*0±|U.(5±). (2.19) 

For aU C e H^°+^{W^) (s > 0), (PH)) and (^1^ allow one to extend ^± as 
a mapping 

VO < s < to + 1, g^ : H'+^/'^{W^) -^ H'-^''^{W^). 

Moreover, the following estimates hold (they correspond to the particular case 
j = of (I23TI1 - (12321) below) 

V0<s<to + 1, |^±7^|ffs-i/2 </i'/*M|<p^|ff., (2.20) 

V0<s<to + l/2, \g^MH^-l/2<^iM\^^p\Hs+ln. (2.21) 

Remark 2.7. From the definition (j2.13p of *p, it follows easily that -i=Q^ re- 
mains uniformly bounded as a first order operator as /i — > 0. However, the 
quantity —Q^ involved in the first equation of (|1.19p remains bounded only if 
we consider ^^ as a second order operator. This fact makes the shallow limit 
singular and requires some specific attention. 

It also follows upon integrating by parts in p.l7p that for all Vi j "02 € 

iJi/2(R'*), 



5± 



V^*0±-P(E±)V''*02±=± / 0i(^±02), (2.22) 
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where (jif and <p2 the (variational) solutions to (|2.17p with Dirichlet condition 
tpi and ip2- This identity is the key ingredient to prove the following estimate 
(see Appendix IB. 1.2p . 

V0<5<to + 1, |(A^^±Vi,A^V2)|<Af^|V'ih.+i/2|V'2|;^.+i/2, (2.23) 

' ' ^^ fj. -"M 

for all V'l) V'2 G H'^~^^^'^{R'^). Finally, we will also need the following commutator 
estimate (see Appendix IB. 1.31 for a proof) 

VO<s<to + l, |([A^^±]7/-l,A^1^2)| <mM|V'i|^-i/2|V'2|^=+i/2, (2.24) 

and the inequality (see Proposition 3.7 of [3]), 

(v-Vu, ^±m) < ^iM\^^\wl,o.\^u\l (2.25) 

that holds for aU v e W^'^^iW^f and u e iJy^(]R'^). 

2.2.3 Invertibility 

It follows quite easily from (|2.22p that Cf^ (resp. t^^) is a symmetric, positive 
(resp. negative) operator. We can also deduce that they are injective (up to 
constants of course) so that the inverse (^^)~^ is well defined on the range of 
g"^ and with values in ii'^/^(]R'^). The next proposition states that {(j^)^^°g^ 
is well defined and uniformly bounded (with respect to e and /i) as a family of 
operators mapping H^ ' {R'^) into itself. 

Proposition 2.8. Let to > d/2 and C e iJ*o+2(]Rd) ^e such that Bl\) is satis- 
fied. For all < s < to + I, the mapping 

is well defined and one has, for all ijj G H^ (^ ), 

where M is as defined in 112. S^) . Moreover, one also has the estimate 

VO < s < to + 1/2, |V(^-)-i o g+ij\Hs < M|V^|ff«. 

Remark 2.9. It is straightforward to check that the proposition still holds if one 
switches the + and — signs everywhere. 

Remark 2.10. Let A : H''+^/'^{R'^) >-^ iJ"-i/2(M'*) be a linear operator such that 

VO < s < to + 1, (AMi/;i, A^i/.2) < mMa(^i)|V'2|^»+i/2, 
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for some positive constant Ma{iPi)- Replacing (|2.23p by this inequality every- 
where in the proof of the Proposition shows that {C}^)^^oA : iJ'^+i/^ -^ ij^^, 
is well defined and 

V0<s<io + 1, \{g-)-^ o Aijj\ ^s+i/2 <MMa{iP). 

For instance, one has, for all 1 < j < rf, 

vo<s<<o + i, \{g-)-^odji^\^s+in < -M|(i + v^|i?|)i/Vlff- 

^ f-i 

Remark 2.11. With a little more work, it is possible to show (using regularizing 
diffeomorphisms in the proof) that for all s > 0, the mapping (g^)^^ o g~^ : 
ijp — > H^ is well defined and that 

\{g-)-' o ^+7M^.+v2 < Mis + i/2)|V'|^.+i/2, 

with M{s + l/2) as in p.3p . but since we do not need such a high order estimate 
here, we do not prove it. 

Proof. Let us prove that ((^^)^^(^+V is well defined in iJ^/^(M''). Wc first prove 
that there exists a unique variational solution $^ £ H^{il^) to the boundary 
value problem 

or equivalently (by Ji2.1.2p . we seek a unique variational solution (j)^ E H^{S^) 
to 

' V^" •F(S-)V^>- =0 in 5", , . 

where E" is an admissible diffeomorphism in the sense of Definition 12.11 and 

satisfying (P?7)) . 

For ah ip£C°°{S-)r\ H^{S-), one deduces from ((2:23)) that 



i^^i^^fU^o) < Mf^m^iM^l^^J^i 



:,\I1. 



Since moreover I'/'i^^olff^/^ ^ A* ^'^||V^ '^\\l'^{s-) (this is a variant of the trace 
lemma, see the proof of Proposition 3.4 of [3]), one has 

The linear form ip i— >■ (^+?/', (^|^^p) is thus continuous on 7J^(5~); moreover, 

V^ •P(E~)V^ is obviously coercive on this space, so that existence/uniqueness 
of a variational solution to (|2.27p follows classically from Lax-Milgram theorem. 
We can thus define "0" — 'Z' |-_o ^*^ *^^* "-"^"^ obviously has g^ip^ — g^'tp and 
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■0 = (^ ) ^ ° Q~^'4'- Since such a solution is obviously unique, the operator 
[g^y^ o g+ is well defined. 

We now turn to prove the estimates given in the statement of the proposition. 
Remarking that \{g~)^^ ° C^+7/'|^s+i/2 = \^]~\f^s+i/2, with V" — 0^| _„ as 
constructed above, we deduce from (I2.19P that 

VJ^Wr^ ° ^+^l^=+V2 < M||A^V''"0-|1l2(5-). (2.28) 

In order to get an estimate on the r.h.s. of this inequality, let us apply A* to 
the l.h.s. of (|2.27p . multiply it by A''^^ and integrate by parts, to get 

/ P{T.-)\7^'' K'cj)- -V"' K'cj)- ^~ ( A%g+tP)A'ij- 
Js- Jr-^ 

Js- 

We can therefore deduce from the coercivity of P(S^), (|2.23p and Cauchy- 
Schwarz inequality that 

/c(S-)||A^V^>-||i.(5_) < MAi|^|^.+v.|^-|^=+i/. 

+ ||[A^P(S]-)]V'^"rllL^(5-)l|A^V^>-||L2(5-). 

Since y^|?/'"|^=!+i/2 < AfJIA^V^ 4>^\\l^{s-) by p.l9p . one can deduce, recalling 
that fc(S-) satisfies ([2J|) . that 

||A^V^>-|U2(5-) < AfVM|7/;|^.+v2 +Af|l[A^P(E-)]V''"riU^(5-)- 
Using the classical commutator estimate 

VO < s < to + 1, |[A^ /]g|2 < |/|ff*o+i|5|ff-i, 
it is easy to deduce from (12.71) that 

||[A^P(E-)]V^>-|U2(5-) < M\\A'-^W^'cb-\\2 (2.29) 

and therefore 

||A^V^"nU2(5-) < Af (VTIIV'l^j+v^ + \\A^-'Wcj,-h). 

The first estimate of the proposition follows therefore from a continuous induc- 
tion on s and (12.281) . The second estimate is then a direct consequence of the 
following lemma. 

Lemma 2.12. Let A be a linear operator mapping H^ into itself for all 

< s < to + 1 and such that \\A\\ rvs+1/2 rvs+1/2 < M . Then one has 

yO<s<to + 1/2, |V(AV')|ff= < Af iV^lff.. 
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Proof of the lemma. From the definition of *p, one deduces that 
Using the assumption made on A, we then get 






< M 



and the result follows from the observation that ,7, _,' .K^a is uniformly bounded 
(as a zero order operator) for fj, S (0, 1). D 

D 

2.2.4 Shape derivatives 

It is known (e.g. [36l|32]) that for all < s < to + 1 and -0± € iJ*+i/2(R''), the 
mapping 

is smooth in a neighborhood of any ( E //*o+2 ^jgd-j satisfying (|2.1I) . Let us denote 
by d^'^±(h)V'=^ (j € N*, h = {hi,...,h,) e if*°+2(K^p) its j-th derivative 
at Q and in the direction h = [hi, . . . ,hj). Such derivatives are called shape 
derivatives. 

A first important result is that there exists an exact formula for the first 
order shape derivative (Theorem 3.20 of [36,), 

dg^{h)i}^ = -eg^{hw^) - £±/i±V • {hV_^), (2.30) 

where 



w 



Remark 2.13. Consistently with the notations used in the introduction, V_ and 
uf^ stand for the horizontal and vertical components of the (nondimensionalized) 
velocity fields in the two layers, evaluated at the interface. 

The formula (I2.30p is not convenient to give control of the operator norm 
of the shape derivatives oi g^ because it looks more singular than expected 
(in fact, the identity sg'^leQw/^ = —e'^/i'^'V ■ Vj^ [HI H] shows that the most 
singular terms cancel one each other). Using direct methods, it is however 
possible to obtain some estimates. Let us first set some notation. 
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Notation 2.14. For all h — {hi, . . . , hj) E X^ {X Banach space), and 1 <l < j, 
we write 

{\y)x^\{Mx and {U)x=\{Mx 

rn=l m^l 

The first estimates are given by (see Proposition 3.3 of [3^ or adapt the proof 
of (|2.33l) given in the Appendix) 

vo < s < io + 1, Wg^{\^)ij\Hs-i,2 < Me'^i^/^{h)H^.t,^l\^\^.^^,2, (2.31) 

and (see Remark 3.3 of ^3, or adapt the proof of (|2.34p given in the Appendix) 

VO < s < to + 1/2, \d^g^ih)i;\H.-i/2 < A/eV(h)H(=+i/2)vto+i|^|^2+i. (2.32) 

For small values of s, (|2.31l) and (I2.32p require much more regularity on the hj 
than on ip. This is the reason why we also need the following two estimates, 
where the s-depending norm in the r.h.s. is on one of the hj rather than on ijj 
(Appendix IB. 2 1 for a proof): for all 1 < / < j, 

\d^g^{h),p\H.-i,2 < AfeV'/^|/Ji|//=+V2(hi)^.v*o+3/2|V|^-*o+i, (2.33) 

for all < s < io + 1/2, and 

\d'g^ {h)ij\Hs-l/2 < M£V|/li|^= + l(hi)^(. + l/2)Vto+3/2|V'|/,(s+l/2)Vto + l, (2.34) 

for all < s < io- 

Remark 2.15. Formulas (|2.31l) and (I2.33P require less regularity on -0 or h than 
([232]) and ((2?34| respectively, but they give a control of size 0{fi^/'^) while ((2?32]) 
and p. 341) give a 0(/i) control. This is reminiscent of the comments made in 
Remark [221 



In the same spirit, is also possible to give two generalizations of (|2.23p to 
shape derivatives (see Appendix IB. 21 for a proof). One has 

|(A^d^'^±(h)^i,A>2)| < MeVl^i|;v.+i/2|V'2|^=+i/2(h)^w*o+i, (2.35) 

for all < s < io + 1, and, for 1 < ? < j, 

|(A^d^^±(h)V'i,A"V2)| < MeV|/ii|H=+i/2(h;)i/.vto+3/2|V'i|^=vto+i|-(/'2|^=+i/2, 

for all < s < io + 1/2. 

Remark 2.16. A consequence of these estimates is that one can use Remark l2.10l 
to prove that (g^)^^ o d^ g^{h) is well defined and give some estimate on its 
operator norm. For instance, one gets with (j2.36p 

\{g-y'^d^ ^ +(h)V'U = + l/2 < e^ Af |/l, l^e+1/2 I V'l /f = Vto + l (h;)^wto+3/2 . 
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2.3 A transmission problem 

Attention is given here to the fohowing transmission problem, whose resolution 
ensures that the velocity potentials in both fluid layers can be recovered from 
the knowledge of ( (i.e. the shape of the interface) and ^. 

(af + ^A)$+ = for - ii+ < z < eC, 

(9^ + /iA)$- = for eC <z<H-, 

5„<i>+|„^, - dn<^-]L, = 0, <9.$±|^^^^^ = 0. 

From the discussion in ^2.1.2[ we can equivalently investigate the following 
straightened version, 

yM± .p(j]±)vM±0± =0 in 5=^, 

^^'^^U-^o-£"<^"^„o=^' (2.36) 

— ^9„o!)^i ^ '^ ~ „ ,-1- 



-dnCf)^\^^^~—-dn(f) 1,^0='^' ^"'^ U = Ti"'^' 



where 9„ denotes as usual the upwards conormal derivative (see (l2.8p V 
Transmission problems of this kind are usually studied with harmonic analysis 
tools and the solution is then given in terms of the single and double layer poten- 
tials (see for instance [H]). Such an approach gives very sharp results in terms 
of regularity but is not adapted to our problem (influence of the boundaries, 
asymptotics, ...). We thus propose an another approach, more elementary, but 
more robust with respect to the applications we have in mind. 

Proposition 2.17. Letto > d/2, < s < <o + l and ( e iJ*°+2(R'^) be such that 
12. 1]) is satisfied. Then for all ip £ i/*"'""'^"(]R'*), there exists a unique solution 
0=^ e i?i(5±) to IfFM) such that A^V*;/)* G ^^(5*) and 



Remark 2.18. As with Proposition 12.81 it is possible to extend the result to all 
s > in Proposition 12 . 1 71 above, and Lemma [2. 191 and Corollaries 12.201 and 12.211 
below. The estimates given in the statement of these results still hold for all 
s > provided that M is replaced by M{s + 1/2) in the right-hand-side. 

Proof. If the result of the proposition holds true, then, denoting ip^ = 'A^U-o' 
one has 

p+V"^ - p"V" = -0 and — -r^+f/'+ = — — ^"V^, 

ri Id 

and thus ^~ — -^(^7 ^)^^^^'0^ (well defined in ij^+i/^ by Proposition 
and 

We now need the following lemma. 
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Lemma 2.19. For all < s < io + I7 the mapping 



^.+ l/2^^d^ ^ iJ^+i/'(R^) 



is one-to-one and onto; moreover, one has 

Proof. Let us first consider the case s — 0. For all ^^ G H/ (K*^), let ip^ = 
^{g^y^ o ^+V'+ e i?y^(M'^) (as provided by Proposition [g^, and let (/)± S 
H^IS"^) be the unique solution of p.l7p . Let us now multiply ^.22^ + by ■jy+P'^^ 
(|2.22p _ by -j^p" and sum up the two identities to get (with P"^ = P(S''=)) 

ii "'5+ IL Js- 

= f ^ig+^+)J[C]i,+ . (2.37) 

Using the coercivity of P^ (with coercivity constants k^ — fc(S*) as in 
(|2.6p l to get a lower bound on the l.h.s., and (12.23^ to derive an upper bound 
for the r.h.s., one obtains 

Now, since V'^ ~ (Q^)^^ ° <? V'", we can deduce from Remark 12.91 that 
\'ijj~^\^i/2 < M|^~Ui/2; using l\2.6\i and p.l9p we then deduce 

1^+1^1/2 <Af|J[C]^+|^i/- 

This tells us that for all < p+ < 1 and p^ = 1 — p+, J[Q is a closed, one-to-one 
operator; it is thus semi-Fredholm. Since moreover J[C] is clearly invertible (by 
a Neumann series expansion) for small enough values of p~ , we deduce from 
the homotopic invariance of the index that it is Fredholm of index zero, and 
thus invertible (since it is one-to-one). We now turn to prove the Lemma for 
the general case < s < to + 1- It is easy to deduce from p.l7p the following 
generalization of ()2.22p . 

Js± 
Proceeding as for the case s = 0, one gets 
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Using (|2.29p and a continuous induction as in the proof of Proposition 12. 8[ we 
then get 

and the resuh foUows from (|2.19l) as in the case s = 0. D 

Thanks to the lemma (and Proposition l2.8p . one can define ■(/'+, -0^ G Hf,, 

by 

V'+ = J[C]-V, i>- = ^{Q-r'og+i^+. (2.38) 

Taking (f)^ as the solution to (|2.17p with -0* as above concludes the proof. D 

As a first corollary to Proposition l2 . 1 7] we can prove that the quantities that 
appear in the equations ()1.19p are well defined. 

Corollary 2.20. Under the assumptions of Proposition \S.n\ and denoting by 
ijj^ the trace of (j>^ at the interface {z = 0}, the mappings 



are well defined and one has, 

VO < s < to + 1/2, |'J^±[£C]V'|h- < M\V^P\h^, 

VO < s < io + 1, Iw^KlV-lff-i/^ < M^/^M|i/'U.+i/2, 
yO<s<to + 1/2, |w±[£C]^lff=-i/2 < nM\ij\^.^+i. 

Proof. This is an immediate consequence of (|2.38p . Lemma [2.191 Proposition 
and Lemma [SHI (and of ([2201) -dMB) for the estimate on w ± [eC]ip). D 



A second important corollary to Proposition 12.171 concerns the operator 
^^[eC] defined in (frT8| . 

Corollary 2.21. Under the assumptions of Proposition \2.17\ the mapping 

^'^ ' ^ ^ g-{p+H+g--p-H-g+y'g+^ 

is well defined for all < s < to + 1/2 and one has 
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Proof. Let us check first that [p'^K^g~ - H_~ P~ Q^)~^ o g+ is well defined. 
We can deduce from the positivity of (p^H^ G~ — p^ H^ Q^) and the identity 

{p+H+g- - p-H-g+) o (^-)-i og+o j[Q-^ = H+g + 

that (p^H^ g^ — p^H_g^)^^ o c^+ is well defined and given by the formula 

{p+H+g- p-H-g+r' o g+ ^ -^{g-)-' o g+ o j[Q-\ 

The estimate given in the statement of the corollary is thus a direct consequence 
of (p:^ and Lemmas Ul^ and Ul% U 

Remark 2.22. Proceeding as in the proof of Corollarv l2 . 2 1 [ one can derive several 
equivalent expressions for (^^[e^], for instance, 

gM] = -^g^-Acr' 

p^gM\ = ■^{i + p-H-g^[eC]oig-)-')g+, 

Remark 2.23. In the undisturbed case (C = 0), it is easy to check by Fourier 
analysis that g"^ = ±(7"! [0, 1] = ±y/JI^\D\ ta.nh{y/p^\D\). It follows therefore 
that 

^ V^|D| tanh(7jF|i^|)tanh(V]r|^ 



Some important features of the operator ^^[e^] are gathered in the following 
proposition. 

Proposition 2.24. Let <o > d/2 and C, e H*"+'^{W^) satisfying iflTTjJ . Then 

(1) The operator -^^[eC] ^^ positive and uniformly coercive on H^ (K'^), 

(2) The bilinear form { — gfj[eQ-,-) is symmetric on H^' x H^' , 

(3) For all Q < s <tQ -\- 1 and i/'i, "02 G H^^ , one has 

\{A'g^[eC]^i,A'^2)\ < AiM|V'i|^.+i/2|V'2|^=+i/ 



1/2. 



Proof Thanks to Remark [121 one has (^^[£C]V'» = ^{g+ip+,J[C]ip+), 
with -0+ = J[C]^^'0. It follows therefore from ()2.37p and (j2.19p that one has 
p\^ip+\ < M{g^[eC]Tl',ip)- Since moreover \'^'ip\ < M\^->jj+\ by Proposition EH 



28 



the first point of the proposition is proved. 

To prove the second point, we proceed as for (|2.37|) to get 



W Js+ Id Js~ 



1^2 , 

with obvious notations for 0^ [j = 1,2). The symmetry of Q follows from the 
symmetry of P* . 

The last point of the proposition follows from the first equivalent expression of 
g^JfC] given in Remark [Z!22l together with ((2?23l) and Lemma [2?T9l D 

3 Symbolic analysis "with tail" of the Dirichlet- 
Neumann operator and consequences 

3.1 Symbolic analysis of {^='' 

It is known that the principal symbol of the Dirichlet- Neumann operator can be 
expressed in terms of the Laplace-Beltrami operator associated to the surface. 
More precisely, \i Q^ [Q denotes the Dirichlet-Neumann on the original (without 
nondimcnsionalization) problem (J1.12p . one has typically 



(3.1) 
where C{s^ C,) depends on the norm of C, in some Sobolev space depending on s. 
We refer for instance to [5TJ |36] for a proof with PDE tools, to j47j for a more 
geometric approach, and to [5] for an interesting paradifferential extension. 
The question that interests us here is the way this identity behaves asymptot- 
ically in the shallow water regime. Since this regime physically corresponds to 
configurations where the effect of the bottom is "felt" at the surface, it is possi- 
ble to guess without computation that p.ip becomes singular in shallow water; 
indeed, the influence of the bottom on ^ + [Q is analytic (by standard elliptic 
theory) and thus not taken into account in the symbolic expansion of (^^[Cjf/; 
(even at the next orders). 

In order to make more precise this singular behavior, let us look at the form 
taken by p.ip in the nondimensionalized setting, and in the case C, — Q. We 
recall that in this case, one has ^■^[O] = ±\//i='=|D| tanh(y7z^|D|), whose prin- 
cipal symbol is -iiyj ^^\D\. It is thus straightforward to check that 



Since ^^[0] and yj ^^\D\ are both of size 0{\/ ^^) (as first order operators) 
while the residual is of size 0(1), the singularity of the symbolic approximation 
of the DN operator is of order 0{ll ^/JJt^). This is the reason why one cannot 
use standard symbolic analysis in a shallow water regime (though this is of 
course legitimate in other regimes, as in [^H01 [T]V 
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Our purpose here is to show that it is however possible to adapt the symbohc 
analysis approach to the present case. As said previously, it is necessary to take 
into account the information coming from the bottom, that is, to include the 
smooth "tail" of the symbol that is neglected in any expansion into homogeneous 
symbols. In the flat case, this "tail" corresponds to (l — ta,nh.{y/ fj,^\D\)) . In 
order to describe the general case, it is necessary to define the usual principal 
symbol g{X, ^) and another symbol t^ {X, ^) associated to its tail; more precisely, 
we define 

5± (X, = ^J^giX, tanh( v/^i± (X, 0), (3.3) 

with 



9{x,o = vW+^M\WW^WoW), 



(x,U - U±e (.)J_^ f + £V(^ + l)'|VCP 

We can then state the main result of this section. 

Theorem 3.1. Let to > d/2 and ( e i/*o+3(R'^) be such that ( fO]) is satisfied. 
Then for allO<s<to and ip G H''+^/^{W'^), one has 

\g^^T S^{x, i?)V'|^.+i/2 < Sfi^^^Mito + 3)|V'|^=+i/2. 

The following inequality also holds 

|^±VT5±(a:,i^)V'|^. < e/xAf (io + 3)|V'|^.+i/2. 

Remark 3.2. In the case d — 1, one has 5(2;,'^) — \^\ and one can compute 

r vi 4±/ c\ ^1 _u ±^Narctan(ey^a^C) 
explicitly P (a;, ^ = 1 ± £ 7^ C- 

Remark 3.3. One readily deduces from the theorem the following estimate for 
the standard symbolic approximation, 

which is the nondimensionalized version of p.ip generalizing p.2p to the nonflat 
case. The interest of the new symbolic expansion of Theorem 13. II is that it gives 
a control of order 0{ey/Ji) of |VV'|_f/s-i/2 (deduced from the 0(£/u'^/'') control 
of I'l/'l As+1/2); it is therefore more precise in its dependence on fj, and e. Even 
outside the shallow water regime (ie if /i = 0(1)), the estimate of the theorem 
is more precise than the standard expansion if e is small. In particular, the 
formula is exact in the flat case (e = 0). Note also that the improvement factor, 
namely Ey/JI = e^^/Ji^, is the typical slope of the wave, a quantity that plays 
an important role in many asymptotic expansions. In the second estimate of 
the theorem, we have a better control of order 0{e^i) but in terms of regularity, 
it is half-a-derivative worse (|^|/^.+i/2 instead of IV'lws)- 
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3.2 Proof of Theorem [XT] 

Since the proof of the result for Q^ follows exactly the same lines as for (^ + , 
we only do it in the latter case. 

We also assume throughout this proof that S(X, z) = {X,z + a{X,z)) is the 
trivial diffeomorphism given in Example 12.31 in particular, 1 + d^a = 1 + e+C 
does not depend on z. 

Let us denote by P = P(X, z,\/,dz) the elliptic operator 

P = V^^ •F(S)V''"'. 
From the explicit expression of i-'(S) given in (|B.7[) . we can decompose P into 

P = P/+£+M+P//, 



with 



^ = i±iiiZ^52_2^+V(T-Va, + /i+(l + 9^fT)A, 



1 + a^cr 

The strategy to prove Theorem 13. II is the following. Since by definition 

0, 



g V^ = a„^|„„, with , ^ _^^ 9„<^|_,-o 

(recall that 9„ stands for the upward conormal derivative, see p.Sp ). we expect 
that at leading order, ^'^ip ~ dn4>app |,^o if 4>app solves the same boundary value 
problem as up to lower order terms. This is the approach used in [3S1 [2] and 
which leads to the standard symbolic analysis of Q'^[Q'<p- The difference here 
is that we want the approximation ipapp to be nonsingular with respect to /i"*". 
This is this extra constraint which imposes the presence of the "tail" in the 
symbolic analysis of the DN operator. 

Since P = P/ up to first order terms that are also of size 0(/i+), the idea is 
to look for an explicit function 4>app satisfying 'Pi4>app — (up to lower order 
terms). In order to construct such a 4)appi we first remark that if we freeze the 
coefficients of P/ and take the Fourier transform with respect to the horizontal 
variables, P/ becomes a second order differential equation with respect to z, 
with characteristic polynomial 

The roots of this polynomial are given by 
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with 



^iX,z,0 = ^ I ^^^^""^p V\e + /^+(|VaP|gP - (VcT ■ 0'), 



l + /i+|Vcr|- 



Kx,z,o = ^^^^v^-e 



We then define (fiapp as the inverse Fourier transform of the solution of this frozen 
coefficient differential operator satisfying the boundary conditions at z = and 
z = — 1, namely, 

(j>app{X,z) = E(X,z,D)ip, 

where the symbol of the pseudodifferential operator S(X, z, D) is given by 

cosh(-\/ /i+ J_j^ a(A, Zj^jdz'j Jz 

The following lemma quantifies the accuracy of the approximation of (p by (j)app- 
Lemma 3.4. The approximation (papp of (j) solves 

Ffpapp = S + fl+h, 

4>app 1,^0 = '0, dn4iapp |,^_i = 0, 

where, for all < s < to, h satisfies 

\\A'h\\2<M{to+3M^s+i/2. 
Proof. One can decompose F(j)app under the form 

P0app = 'Pli'app + e'^t^'^^Ili'app] (3.4) 

we now analyze the two components of the right-hand-side separately. 
- Analysis ofFjc/fapp. From the relations 

Va-Vd,<l)app - Vcr • VOp(a^S)V' 

{l + d,a)A^app - -0p((l + a,a)5|enV' + (l + 5.a)[A,0p(S)]^, 
we deduce that 

P/</>app - Op(P/(X,z,ie,a,)S)V'-2A*Op(Va-V9,S)V'+A*(l+9.a)[A,Op(S)]V'. 
From the definition of S, it is also easy to check that 

2 /-O 



P,(X,z,»e,a.)S = e+^+i±i^^^!expHM+/ &) 



a^a sinh(y^/^^a) 9^5 cosh(y^/^^ a) 

+ i 



/i+ cosh(-\//i+ J -1^ a) £ cosh(-\//i+ J j^ a) 
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since moreover 

(1 + 9,a)[A, Op(S)] - Op((l + d,<j)AE)) + 20p(i(l + a,a)VS • 0, 
one easily gets that 

Pi<Papp = e+/i+Op(i?i + i?2)V' + £+M+(l + a,a)Op(^)V', (3.5) 



with 



2 /-o 



1 + o^a J^ 



id.h 



' IJL^ cosh(-\//i+ J-^ a) cosh(-\//i+ /-^ a) 

"i?2(^,2,0 = 2i(l + a^cr)VS-C-2Vcr- V9^S. 



We now turn to control all the components of p.Sp : 
- Control of Op{Rj)t/j, (j ~ l,2j. Let us remark that 

Op(i?,)^ = Op(R,) • V^t ^j ^ 1^2, 3), 

with tp^ = exp(co-\//i+z|I?|)'i/; and 

R,(X,z,e) = -^^^^^^^||exp(-coVM+z|e|), 

and where cq > is such that i^f^^iy^ii > 2co for all {X, z) £ S (and thus 

a(X, z,^) > 2co|'^|). It follows from the definition of Rj and cq that Hj is a 
pseudo-differential operator of order zero, whose coefficients depends on second 
and lower order derivatives of ( (through a). Using Theorem 1 of [37], which 
gives precise estimates of the operator norm of pseudo-differential estimates, we 
then get that 

VO<s<to, Vze[-1,0], \Rj{X,z,D)f\H^<M\f\H^ 

(it is easy to check that the possibly singular terms l/e^ in the definition of Rj 
are compensated by the e+ contained in ct = e+(z + I)(^). We thus deduce that 



|A'''Op(i?i+i?2)V'l|2 < Af||A"V^ 



2, 



< 



M|^|^.+i/2, (3.6) 



the last inequality following from an easy computation (recall that z < 0), as in 

Proposition 2.2 of [3] for instance. 

- One can prove along the same lines that 

||A^Op((H-a,a)^)^||2<Af(io + 3)|VL.+i/2 (3.7) 

(there is Af (to + 3) rather than M in the r.h.s. because AS involves third order 
derivatives of Q . 
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From (133]), dSUl) and dXT]), we get that 

Pi(f>app -= e+n+hi with WA'hih < M{to + 3)|V^|^.+i/2. (3.8) 

- Analysis of'Pii(j}app. From the definition of P// and (j)app, one gets easily 

P//0app = Op(i?3)V', with i?3(^,^,C) - ^{—^d,\Va\^ - Aa)9,S. 
Proceeding as in the previous point, we get, for all < s < to: 

||P/7<^app||2<M|^|^. + l/2. (3.9) 

- End of the proof. The result follows directly from (IXi)) . (P?5| and (jX^ . with 
h = hi +Pii(t)app (it is straightforward to check that 4>app satisfies the boundary 
conditions) . D 



(3.10) 



It follows from the lemma that the difference u = (p — <f>app solves 
yM+ . p(E±)V''^M = -e+fi+h 

The following lemma allows some control on u. 
Lemma 3.5. Let u e H^{S) and h e L'^{S) be such that 

V^ •F(S+)V^^u = /i, 

Then, for all < s < to, 

||A^+iV^^u||2 < -^M\\A'h\\2. 

Proof. Let us first remark that ft, = V^ • g, with g = {J^^ h^e^ and we thus 
deduce from Lemma [B. II that 

||A^V^\|l2 < M|lA«g|l2 < M\\A'h\\2. (3.11) 

For 1 < J < rf, one remarks further that v = djU solves 

V^^ • P(S+)V^^t; = V"^ • g, 

with g = -djP{T.+ )V^u + -^JiBj. We thus deduce from ((RT|) that 



||A'^9jV''"u||2 < M||A^g||2 

< M(l + ^)|jA^ft||, (3.12) 

where we used (|3.1ip and (jB.2p for the second inequality. The lemma is thus 
proved. D 



34 



Applying Lemma [531 to p.lOp with h = — e+^+/i, we get (replacing e+ and 
/x+ by e and /i since the ratios H^ / H^ and H^ /H~^ arc controlled by M), 

< ey^M{to + 3M^.+i/2, (3.13) 



where the last inequality follows from Lemma 13.41 

In order to prove the theorem, let us now remark that 

By Green's formula, we get therefore, for all (p G L^(M'^), 

= / P(E+)V^\ • V^+A^+i/V^ - e+fi+ f hK'+^/\^ 
Js JS 

= [ A^+ip(S]+)V'^^u • A-i/2vA.+ <^t _ ^+^+ f A^/^Al/V^ 

J5 Js 

where (p'^ is defined as 

(x being a smooth compactly supported function equal to 1 in a neighborhood of 
the origin). Since \\A-^/^Wt''' ip^\2 < iiJ+)^^'^\(p\2 and ||Ai/2^t||2 < (^+)-i/4|<^|2 
(see Prop. 2.2 of [3]), we get therefore 

A^+i/2(g+^ _ a„^„^^ i^^J^ < M{ti'^^\\A'+'V^\\\2 + e^i'/^A'hh)\ip\2 

< £/i3/4Af(io+3)|VU. + i/2|(^|2, 

the last line being a consequence of p.l3p and Lemma 13.41 Since moreover one 
computes easily that 



dn<f>appU^„ = e, •P(E+)V^\ 



'app l^^o 



M+VleP + M+(|VCP|CP - (VC • 0') tanh( / a). 



the first estimate of the theorem follows from a standard duality argument. 
For the second estimate, we proceed as above to get 

/ A^(^+V'-a„0appu_o)'^= /A^+ip(S+)V^^u-A-iV^^^t_£+^,+ /'a 



'hifK 



Since ||A -^V^ (p^||2 < y^|(^|2 and ||(y3^||2 ^ Iv'b, we deduce as above by a 
duality argument that 

|^+V-9„0appU=olH'- < V^A/f|lA^+iV^''u|l2 + e/i|lA^/i|l2 

< e^M(to + 3)|?/'Us+i/2, 
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where we used Lemma [33] (with h = — £+/i+/i) and Lemma [?^ to get the second 
inequahty. The end of the proof is then exactly as for the first estimate. 

3.3 Symbolic analysis of (g^y^g^ 

Theorem 13.11 shows that Op(5'''), with 5* as given by p. 31) . provides a good 
description oi (^^ . It is thus natural to expect that a good symbolic description 
of the operator (^^)^^^^ (which exists by Proposition 12. 8p is given by 



(^-)-i^+^Op(— ) = 0p( 



5'+, ^ / H+ tanh{y/Jj+t+) 



H tanh(-\//i^i^) 



with i* (X, ^) as in Theorem 13.11 The following corollary shows that this is 
indeed the case. 

Corollary 3.6. Let to > d/2 and C G H*°+^{R'^) be such that ![2l\) is satisfied. 
Then for allO < s <to and ip G iJ^^^/^^gd^^ ^^^ ^^^^ ^-^j^ k^O,l, 

\[{g-)-'g+-Op{^)]i;\^.^+U2<efi-''/^Mito + 3MHs-./2. 

Remark 3.7. As for the symbolic analysis of (J' ^, it is crucial to take into account 
the nonhomogeneous tail of the symbol. A simple look at the flat case shows 
indeed that one cannot expect an estimate better than 0{l/fj,^^'^) if the principal 
symbol (constant equal to —H^/H~ here) is used instead of S' + /S~ . In the 
case fc = 1, e have the same 0{ey/JI) improvement as in Theorem 13.11 Note 
also that the /i^^/"* in the r.h.s. of the estimate given in the corollary is not 
singular since it corresponds to the fact that the H^ -norm (rather than the 
TJ'^+^/^-norm) is considered in the l.h.s. In the case fc = 0, there is a further 
gain of fj,^^'^ in the estimate, but a loss of half-a-derivative, as in Theorem 13. II 
Remark 3.8. Contrary to Proposition 12.81 for instance, Corollary 13.61 requires 
that tp belongs to a standard Sobolev space rather than a homogeneous one. 

Proof. Denoting by S^ {X, ^) the symbol of the DN operator (^ * as given by 
Theorem 13. 1[ we define 'tp~ = Op (4^ )-(/'. We also denote by 0~ the solution of 
(|2.17P - with Dirichlet condition V'~ at z = 0. Then, by definition oi g^, one 
has dn4>7 _ = Q~4'~ ■ If 0~ is as in ()2.27p . the quantity u = (j)~ — (p^ solves 
therefore 

I 5„M|^^„ = g^^p^ - g~^f~, dnu\z=~i = 0. 

Proceeding as in the proof of Proposition l2.81 we get that 

/ P(E-)V^"A"u- V"A"u = - I K''{g+'tp~ g'i}^)K'uQ 

+ [ [A^F(S-)]V''"w•V^"A"u, 
Js- 

where uq stands for u\^^g. Let us now state the following lemma. 
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Lemma 3.9. One has, with k = 0,1, 

Proof of the lemma. We first consider the case k = 1. Let us remarllj that 

with 

Ri = ([A^^+]VAA«uo) 

i?2 = (A«-i/V,Ai/2(^+-Op(5+))A«uo). 

Remarking that ^^/'*|*PVI//=-i ^ IV'Iff^-i/a, we can use (|2.24l) and Theorem 13. II 
to get 

l^jl < liM—^\lP\Hs-l/2\uo\j^s^+l/2 (j = 1,2). 

RecaUing that i/i^ = Op{S^ /S^)4>, we write similarly 

(A^(^-V^-,A^z.o) = (Op(|j)*Op(|^)A^V,A^'p2^o)+i?'i + i?^, 
with 

i?'i - ([A^^-]V^-,A^uo) 

i?^ = (A^-i/2v;-,Ai/2(^-_Op(5-))A^uo), 

i?:, = (AV2[A^Op(|^)]V^,A-V2op(|_)A^^2^„). 

Remarking (this is a consequence of Theorem 1 and Corollary 30 of [37]) that 
||0p(iS'+/5^)||^s-i/2^/^s-i/2 < M {0 < s < to), we get as above that R'l and 
i?2 satisfy the same estimates as Ri and i?2 respectively. Using Theorems 1 and 
8 of [37], we also get that 

\[A^,0p{S+/S-MhU2 < eM\^\Hs-^;2 

since /Li|*p^uo|//s-i/2 < M^^^I"o| r>^+i/2, this allows us to conclude that R'^ satisfies 
the same estimate as Rj and i?' (j = 1, 2). 
We thus have 

(A^(^+V - g'i^-),A-'uo) =Ri + R2 + R[ + R'2 + R3 
+ (A^/'(Op(|J)*-Op(|^)Op(|^))A>,A-V2^2^^) (3.14) 



^We use the same notation 53 for the operator j—^ — r-rrr and its symbol 15 , . ,„ 
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We then deduce from the above estimates on the residual terms that 



I|0P(^)* - 0P(^)0P( — )IIh-V2^,,i/2 < e/iM, 



The result follows therefore from the observation that 

which is a consequence of the composition estimatqil of Theorem 8 of 13^ and 

of the estimates |10p(|i)* - Op{^)\\h-i/2^h^/^ < ^M^^- 
We now briefly indicate the modifications to be performed in the case fc = 0. 
For the estimate of Ri, we just replace the control ^^/^|*PVIh=-i ^ IV'lff»-i/2 
by I^VIff'-i ^ IV'Iff"- For R2, we rather write 

i?2-(A>,(^+-Op(^+))A^wo), 

and use Cauchy-Schwarz inequality and the second estimate of Theorem 13.11 
The same modifications must be done for R[ and R'2 respectively. For R'^, we 
write 

R', = (A[A^Op(|^)]^,A-lOp(|^)A^^2^o). 
and use the estimates 

\[A^Op{S+/S-Mh^ < eM\i,\H^ 

The control of the last term in (|3.14l) is modified along the same lines as R'^ and 
the result follows. D 



Thanks to the lemma, one can proceed as in the proof of Proposition 12.81 to 
get (in the case k — 1, the modifications for the case k = are straightforward), 

IIA^V^^Ib < eM{to + 3)(m^/4|V'|^.-i/2 + WA'-'V^"' c^'h), 
and thus, after a continuous induction on s, 

IIA^V^^Ib < e/i'/'M(to + 3)1^1^.-1/2. 
The result follows therefore from (|2.19p and the fact that uq — ((§'^)^^(7* — 

Opis+/s-))ip. n 



^Theorem 8 of |37) deals with commutator rather than composition estimates. However, 
the commutator estimates of Theorem 8 follows from a composition estimate, exactly in the 
same way as Theorem 7(iii) follows from Theorem 7(1). 

**This estimate on the adjoint is not stated in |36| . However it can classically be derived 
with the same techniques as the commutator estimates of Theorem 8 of that reference. See 
also Proposition 1.8 of |26) and Chapter 13, §9 of [51) . 
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3.4 Symbolic analysis of {g-)-^g^[eC] 

The following corollary shows that the asymptotic description of {g ~)~"^^^[eC] 
that one naturally expects from the definition of (^^[eC] in Corollary 12.211 and 
Theorem 13.11 is correct. In order to state the corollary, it is convenient to 
introduce the zero-th order symbol 

S,iX.O^,*-,-^W^ (3,15) 

Corollary 3.10. Let to > d/2 and C G i/*"+3(R'') be such that l2l\) is satisfied. 
Then for all < s < to and %}) G H''~'^/'^{M.'^), one has, with fc == 0, 1, 

|[(^-r'^MK]-;^Op(^)]^|^.+v. <£M-'/'M(io + 3)|Vli,=-./=. 
Proof. Thanks to Remark I2.22[ we can write 

It is a direct consequence of Corollary 13.61 that Op(S'j) provides a good approx- 
imation of J[Q. The following lemma shows a corresponding property for its 
inverse. Throughout this proof, we only consider the case fc = 1; the modifica- 
tions that must be performed in the case fc = are absolutely similar to those 
made in the proof of Corollary l3.6l 



Lemma 3.11. Under the assumptions of the corollary, one has, for all < s < 
to and fc = 0, 1, 

\{J[C]-'-Opi^)iJ\^.+^;.<e^l-''^^Mito+3MH^-.r2. 

Proof of the lemma. From Lemma 12.191 we know that 

|(J[C]-^ - Op(^)7^|^.+i/. < M| (1 - J[C]Op(^))V|^.+i/. 



< 



Af (1(1 - 0p(5,/)0p(^))^|^.+v. + |(J[C] - 0p(5j))0p(^))^|^.+v2) 
M{to + 3)(47ll(l - Op(^j)Op(^))V'|^.+v. + ^|Op(^)V^lff=-i/0' 



where we used Corollary |3.6l to derive the last inequality. Since the pseudodiffer- 
ential estimates of '36 show that ||1 — Op(S'j)Op(l/5'j)||^s-i/2_^j:^s+i/2 < eM 
and ||Op(l/5,/)||^s-i/2^//s-i/2 < M, the result follows easily. D 

To conclude the proof of the corollary, let us remark that | [{g^)^^gfj.[£C] ~ 
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-^Op( g_g )]V'|jj^+i/2 is bounded from above by 

< ^\[ig-)-'g+Opij-) - op(^^)] vl^.+i/. + j^M{to + 3Mhs~u2, 

where the second inequahty comes from Proposition l2.8l and Lemma [3. Ill It is 
then an easy consequence of Corollarv l3.6l and of the pseudodifferential estimates 
used many times in this section that the first component in the r.h.s. of the last 
equality is also bounded from above by ^74 Af (io + 3)|'0|ifs-i/2, which concludes 
the proof. D 

3.5 Symbolic analysis of {p^-^g^ — p^ ^Q^Y^dj 



For the sake of clarity, we write 

Q=-iY^g--E'^Q+) and S=-{p+^S--r^S+). 

~ ^ ~ ~ ~ (3.16) 

The operator "^"^oQ -^dj (1 < j < d) is well defined thanks Remark [2TT0I It will 
be useful to describe it using a pseudodifferential operator, as in the following 
corollary. 

Corollary 3.12. Letto > d/2 and C e i/*"+3(M'^) be such that |gj]) is satisfied. 
Then for allO<s<to,k = 0, 1, and f £ H'+^^'^{R'^), one has 

\V' o g-^d,f - Op(£)9,/|^.+,/, < ^i^Af(io + 3)1(1 + ^\D\)'^'f\^^. 
Proof. We can write 

'P^ o g-^d,f = Op(^)a,/ + R1 + R2 



where, using the notation f ~ Q ^djf 



1 Op(^) o op(A)l^V, R2 = Op(^) [Op{S) -g]f. 



s ' ^*P- 

From the composition estimates of 13 7j . we get, with A: = 0, 1, 

|i?i|ff=+M2 <eAf(io + 3)|<p2/|^_,^,,, <e/i~''^^M(io + 3)|W|ff=, 

the last inequality following from the observation that 1*^/1^:^1/2 < M~^^''l/lffi 
and |^/|2 ^ l/lffi- Since moreover Remark 12.101 implies that 

IWk=<Ai-'M|(l + VA^|i?|)i/2/|H=, 
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we deduce that Ri satisfies the bound given in the statement of the corollary. 
For i?2, we get 

\R2\h'+^/2 < -M{to + 3)\[OpiS)~g]f\Hs+.,2 

the second inequality being a consequence of the first {k = 1) and second {k = 0) 
points of Theorem 13.11 One then shows as for Ri that i?2 satisfies the desired 
estimate. D 



4 Quasilinearization of the equations 

The goal of this section is to derive a quasilinear system from the two-fluid 
equations (J1.19p . Throughout this section, the operators §'^[eC] and (^^± [e^C; 1] 
are denoted (^ and (^ ^ respectively, when no confusion is possible. 
For all j eW, we also denote by d^ g {h)4> and d^g^(h)ip (h = (/ii, . . . , hj)) 



respectively in the direction h. 



the j-th order derivative of the mappings ( i-)- g^leQ'ip and C '~^ ^ ± [s Cj 1]V' 



4.1 Definition and basic properties of some new operators 

We define first a linear operator T [U] of order one, 

TPU = V • ifv+) + p-H-g o ig-)-'{v ■ (/K^D), (4.i) 

where V_ and ur^ are defined using the mappings w * [eC] and 1/^ [e(^] of Corol- 
lary [2^ as 

y± ^ n/^lsC]-!!! ~ ew±VC, w± = fv^[eC]ip. (4.2) 

In the water waves case {p~ = 0), this operator simplify considerably into 
"^[Ujf = '^ ■ (/I£^)- The symbolic analysis of Section |3] allows us to prove that, 
in the general case, it behaves roughly the same way in the sense that it is a first 
order operator with skew symmetric principal symbol. Note that it is necessary 
to use here a symbolic analysis "with tail" as in the previous section in order 
to have a non singular dependence on the parameters in the estimates of the 
proposition. 

Proposition 4.1. Let to > d/2 and U = {C,ip), with ( e H^°+^{R'') satisfying 
fmp and iP e iJ*«+3(K''). Then 

(1) For allO<s<ta and f e H''+^/'^{R'^), one has 

(2) For alla=l + b, with b e H^«+^{R'^), and all f e L^{R'^), 

(aT [U]f, f) < M{to + 3)(1 + |b|H'o+0|VVU*o+i l/li 
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(3) Let K be adxd matrix with entries in iJ*o+i(R''); then for all f € H'^{R'^), 

(V • KVT [[/]/, /) < M{to + 3)|XU*o+i |VV|i,*o+i I/I'hi . 

(4) For all f e iyi/2(R'i), g e H^/^{R'i), one has 

(^[f/]*/,3)<A^|VV|ff*o+2|W|2(|5|2+/i'/'|<?|ffi/2). 

Proof. In the proofs below, we reduce 1 \U] to its second component in \4-l% 
since the estimates are much easier to prove for the first component and can 
therefore be omitted. 
For the first estimate, we write Q o {g^y^V ■ {flV_^l) = A + B with 

,,i/4|r)|i/2 

A = goig-r^v-i^f-^j^ifiv^)) 

Recalling that g = -i-^g'+ o J[(]~\ we use (fZ!20l) . Lemma [2?T9l and Remark 

12. 101 to get |A|^a-i/2 < A/|/jy''']|^s+i/2. Proceeding in the same way, but using 
(|2.2ip instead of (I2.20p . one can chech that B can be controled similarly. Since 
moreover, with < s < io, one has |/|Vl''']|/fs+i/2 < Af |/|^s+i/2|VV'|//to+i/2 
(from the definition of]/* and Corollary 12. 20p . the result follows. 
For the second point of the proposition, we write 



(ac^ o {g-)-'^ ■ iflV^j), /) = (V • (/|y±]),Op(^)(a/)) 

+{flv%vKg-)-' o g op(^)](a/)). 

Since the operator / i-> aOp( g'^g )*V-(/|Vl^|) is skew symmetric at leading or- 
der, the result for the first term of the r.h.s. follows from the pseudo-differential 
estimates of [371. In order to control the second term, we decompose it into 
Ci -I- C2 with 

,,l/4|r)|l/2v7 c+ 

^1 = (/K^i i + j/4|^|i/2 ((gT^°g-Qp(^^))(^/))^ 
C, = {fiv% ^^ ^J^^^^^^^ iig-)-' o g -Opi^maf)). 

From Cauchy-Schwarz inequality and Corollary 13.101 (with fc = 1), we get 
|Ci| < M{to + 3)|/|[V:'^l|2|o/|2. For C2, we proceed similarly but with fc = in 
Corollary 13. 101 to get that C2 can be controled as Ci and the result follows. 
The proof of the third point being yery similar to the second one, we omit it. 
For the last point, the term to control is 
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It is therefore a consequence of Cauchy-Schwarz inequality, Proposition 12.81 and 
Lemma [231 that 1^1 < Af IWbKl + y^lDD^^^iglV'^DU, from which one de- 
duces the resuh by a standard commutator estimate and Corollary 12.201 D 

The linearization formula for d"' {g ^[eQip) (see Proposition 14.31 below) is 
quite simple if given in terms of the "good unknowns" C(a) a-nd V'(a): where ('(•q,-) 
and -0(0) are defined as 

C(a) = 9"C, (4.3) 

V'(a) — d'^ip — ewd°'(, with w — p'^w'^ — p^w^, (4.4) 

or equivalently, if '0* = ^i _ (0^ solving (|2.36p ). 

The principal part of the linearization formula for 8°" {C^ fj_[e(]tp) is given in terms 
of Q^(a) and T \U]d"'C,, but to handle the surface tension, it is also necessary to 
take into account the subprincipal part (with respect to ip). In order to describe 
this subprincipal part, it is convenient to introduce the following notation, 

C(q> = (C(q1),- ■ • ^Cca-i+i)), "0(0) = ("(/"(qI), ■ • ■ , V'(q<J+i)), (4-5) 

where a^ G N'*+^ is such that a^ + ajGj = a. 

We can now define the operator Q{a)i^{a) that arises in the description of the 

subprincipal part of d°'{Q^[eQ'\jj), 

d+l 

Q(a)ip{&) =^ajdg{dj0^p^aj). (4.6) 

Next, the linearization of the surface tension term shows that 

where the dots stand for derivatives of ( of order lower or equal to jaj and C(c) 
is as in (|4.5p . The positive definite d x d matrix JC [VC] is defined as 

_ (i + |vcp)id -VC0VC f, 7. 

while the second order operator ^(q.) [VC] is given for all i^ = (/i, . . . , fd+i)"^ by 

d+l 

^(a)[VC]F = -V • [^ {d^{Wd,C)S/f, + dJC{Vf,)Vd,<:)]. (4.8) 

Finally, we introduce here the first order operator £ [C] : iJ^/^(R'*)'^ — >■ 
iJ^^''^(M'*)'^ that plays an important role for Kelvin-Helmholtz instabilities, 

'E[C] = Voip+^g--p-^g+)-^oV^; (4.9) 

the following proposition gathers some of its main properties. 
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Proposition 4.2. Let to > d/2 and C € H*"+'^(R'^) be such that l[KT]) is satis- 
fied. 

(1) There exists a constant c< M such that for all V G iJ^/^(R''), 

Q<^,['E[(:]V,V)<c\[l + ^\D\Y'^V\l■ 

we denote by e(C) the smallest such constant. 

(2) If moreover C, is time dependent and satisfy \2.1\l uniformly in time, then, 
for all F e H^/'^{W^Y 

\{[du'E[C]]F,F)\<'-M\d,CU\{l + ^\D\Y'^F\\. 
M 

Proof. With Q as in p.l6p . we can write (£[C]T^, V^) — [Qw,w), with w = 
Q~^V -V and the positiveness of £ [^] foUows therefore from the positiveness of 
Q . Since moreover Q — —-^Q^ ° J[C\^ the upper bound fohows from Lemma 
EH and Remark [Sini ~ 

For the second point of the proof, we remark that [dt, Q^^] = ^§^^[dt, Q]Q~'^ , 
so that _ 

{[du^[C,]]F,F) = -{[dug]g~^v-F,g-^v-F). 

Using Proposition 3.6 of j3], we deduce that 

\{[dt,'E[C\]F,F)\< ^l^mg-^v ■ F\2, 

and the result follows from Remark l2.10l □ 

4.2 Linearization formulas 

We establish here linearization formulas for the operator g . These formula play 
an important role in the quasilinearization of the equations presented in Section 
14.31 Proposition 14.31 below gives a linearization formula at leading order; in 
order to measure the size of the neglected terms, it is necessary to introduce at 
this point the energy S^{U) defined for all iV G N, cr > 0, and U = {(, ip) as 

£j^(t/) = |VV.|^.„+. + Y. l^"CI?,i+IV^(a)l|y.- (4.10) 

Q,gN''+Ma|<Af 

It is convenient to introduce a specific notation for the terms involving C, only 
in this energy, namely, 

icIh-«= E i^"ci//i (4.11) 

aeN''+i,|Q|<Ar 

(the difference between H^'^ and H^'^^ is that the former allows for time 
derivatives of C)- Another quantity that will be related to the energy through 
the stability criterion derived in Section [5] is |CI<Af+i/2> defined as 

ICI<iv+i/2>= E W^d%\^. (4.12) 
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For all T > 0, we also denote by -E^ t ^^e functional space 

E^T = {Ue Ci[0,T];H^ x H*°+^R'')), sup £^{U{t)) < oo}. (4.13) 

o<t<r 

We finally denote by m^{U) any constant of the form 

m^{U)^C{M,£^{U),:^), (4.14) 

with M as in (j2.2p . We can now give an important linearization formula for 
g (note that for d = 1,2, the condition (|4.15l) is satisfied with to = 3/2 and 

N>5). 

Proposition 4.3. Let T > 0, to > d/2 and N eN be such that 

[(A^ + l)/2] > 1 Vto + 1/2 and N>tn + 7/2. (4.15) 

Let also U = (C, ip) E E^rp he such that 12. 1\] is uniformly satisfied on [0, T] . 
Then for all a e N''+^ with 1 < |a| < N, one has 

-a"(^7A) = -^V(«)+ei?a, (|a|<iV-l) 
^ /i 

-9"(^7A) = -gi^io.) - eT [UWC + -^(a)^(d) + £i?a, (|a| = N) 

where the linear operators "T [U] and Q (a) o-f£ defined in |^. j[ j and |/^.5[ j while 
Ra satisfies the estimate 

yO<t<T, \R^{t)\Hi<m^{U{t)). 

Proof. Let us denote by (p^ the solution of the transmission problem p.36p 

provided by Proposition 12.171 ^nd let us write ip^ — 1 ^^ (or equivalently 

define ip^ as in (|2.38p ). For the sake of notational simplicity, we also omit 
to write the dependence on the time variable when this does not raise any 
confusion. Consistently with the notations of §2.2.41 '^e also denote by dg {h)ip 
the derivative of the mapping ^ i— >■ CJ'^[eC]'0 in the direction h. 
We focus on the most difficult case, namely |a| = A^ (when |a| < N — 1, it is a 
consequence of the first estimate of Proposition 14. II that e'T [U]d"( can be put 
in the residual term). The proof relies heavily on the following lemma (note 
that the last estimate of the third point is not used in the present proof but will 
be necessary to establish Proposition 14. 5p . 



Lemma 4.4. Let to > d/2 and C € H*'>+'^{W'-) be such that l[Kl]} is satisfied. 

(1) For all ij) £ H^/^{R'^) and h e iJ*"+2(M'*), one has 

dg {h)i; = -eg (hlp^ur^i) - e/iT [U]h. 

(2) One can replace g^ by g or w=^[eC] in ^E3J\) - [£34^ . 
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(3) Let r > and S e N''+^ be such that r + \S\ < N - 1, with N as in the 
Proposition. Then one has 

iqja^Vlffj+i < m^(C/) and \^d'^\Hr+i < m"" {U){l + efi^^^\CUN+i/2>)- 

Proof of the lemma. (1) Recalling that § — jj-f'j^ ° J[C]~^j where J[(] is as in 
Lemma [533 one gets 

dgih) = -l^dg+ih)j[C]-' - g o dJ[C]ih) o j[c]-\ 
ti 

Since moreover one can compute 

dJ[(:]{h)^^-^{g-)-'odg-{h)o{g-)-^og+-^-^{g-)-^odg+{h), 

we get, observing that i/'^ — >^[C] "^V' and ?/'^ = ^ffriQ)^^ ° g^i'^i that 

dg{h)i, = 2^dg+{h)^+-p-g o ig-)-'{dg-ih)^- - ^dg+{h)^+). 

We can thus use ()2.30p to get 

dg{h)^ = %-{l + p^ H~ g o {g^r^)g^{hw+) + ep- g (hw-) - efiT[U]h. 

H 

The result is then a consequence of Remark 12.221 

(2) Proceeding as in the first point, it is possible to show that d^ g (h) and 
dw^{h.) in terms of C^ + , o/^[C]~^ and their shape derivatives, so that one can 
deduce the result from ((23T|) - (|234)) . Remark [2l6l and Lemma [2l9l 

(3) Let us first prove the estimate on the iJ^+^-norm of ^d^tp. Thanks to the 
assumptions made on r and S, one has 

;3eN''+i,|^|<Ar-l 

Since moreover d"ip = ip(j3) + ewd^C, we get 

,3eN'i+M^|<jv-i 

< max{l,— =} Y^ (|<PV(a)|2) +£Nffto+2|ClHi."), 

^ QeN<i+MQ|<JV 

where we used the identity |^/|2 < l^/b to derive the second inequality. Since 
w = p'^w~^[C]ip — p~w~[C]tp, the result follows easily from Corollarv 12.201 
For the estimate in iy^^-norm, one derives along the same lines the identity 

|q39Vk-+i< E (IW(«)|2 + £|q3M"C)|2). 

aeN''+'^,\a\ = N 
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From the product estimate (that fohows directly from Lemma 5.1 of [3? for 
instance) , 

l/ffl^v^ < |V/|h.o \g\2 + |/|oo|5Uv^, (4.16) 



and the observation that |w|oo < v^*^ (^) (Corollarv l2.20p . we get I^Pd^c^^Ob 
m^(C/)(l + ey7Z|<pa"C|2). The resuh then fohows easily. D 

Denoting by {9", C^IV^ the Poisson bracket (with a^ as in ()4.6p ) 

d+\ 

we can write 

d'^g^ = ^9"V + rf^ (a"C)^ + ^(a)^(a> + e^^a, (4.17) 

with EfxRa = ({9", ^ IV' — Q{a)i^{a) ) + SM-^a ^md where e^i?^ is a sum of terms 
of the form 

d^g{d'\,...,d'\)d'^=:A,^s,., (4.18) 

where j e N, i = (i\ . . . , t^) G N^'^+i'j' and 5 G N'^+i satisfy 

i 

^|i*| + |,5| =iV, 0<|(5|<7V-2 and \i'\<N. (4.19) 

We first give some control on the Aj^s,^ (and thus on R'^)- Let / be such that 
|i'| = maxi<j„<j |(.™|. Let us distinguish four cases: 

(i) The case \5\ = j = 1 and l — i^ (and thus \l\ — N — 1). One then has 
^j,s,L = dg{d''C)d^tlj, and we can therefore use the first point of Lemma l¥^ to 



< 



get 






We can thus deduce from Corollaries 12.201 an d2.21[ and Proposition 14. 1 1 that 
-\Aj^5Ahi < eM(to + 7/2)|C|^i,«|*P9'v^|^*o+V2; 

since io + 7/2 < N, these two estimates, together with the third point of Lemma 
[441 yield that 

\Aj,sAhi < efim'^iU). (4.20) 

(a) The case \5\ < N -IVto- 3/2 and |t'| < iV - 1. It follows from ([2341^1 
that 

-\A,,s,]Hi<eM\d''C\Hy- n |5''"CI^-*o+5/2|*P9V^|^lV*„+3/2. 



^We still refer to 112.3411 if (J^ is replaced by §, thanks to Lemma 14. 4l ii: we do the same 
for all the identities involved in Lemma 14. 4l ii 
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From the assumption that |i'| < TV — 1 and the definition of Z, we also get that 
for all m y^ I, \l"^\ < [{N — 2)/2], and thus, from the assumption made on TV, 
we have |t™| + 1 V to + 3/2 < N. Using the assumption on S, we deduce from 
the last point of Lemma 14.41 that (|4.20p holds. 

(in) The case \5\ = and |t'| = A^ — 1 (we take for simplicity I = 1). Then 
one has Aj^s,i, = (fi§{d'' (, 9' Qip. Inequality (12.33^ provides a good control in 
terms of regularity but not with respect to its dependence on fi. We can remark 
however, after differentiating the formula stated in Lemma U^i that 

-ey/JIg ihilp^^dw^{h2)^) - e^idl {h2)hi. 

With hj — d'' (, one can use (|2.3ip to control the first term of the r.h.s.. Corol- 
lary 12.211 and Lemma I4.4l ii for the second one, and Lemma I4.4l ii for the last 
one. Inequality (|4.20l) is then satisfied. 

(iv) The case \5\ > N -- I V to — 1/2. Using (|2.32p and the second point of the 
Lemma, we now get 

j 

m—l 

Since, by (I4.19p . |(^| + 2 < iV and \l"^\ < io + 2, one can conclude as in the 
previous case that (|4.20p holds for iV as in the statement of the proposition. 

We have therefore proved that the residual R'^ satisfies the estimate stated 
in the proposition. Remark now that -{{d", g}ip— ^(q)V'(q)) is a sum of terms 
of the form (with w as in ()4.4p ) 

-^dgidX)i^wd\) (H = i, L + 6^a) 

that are bounded in _ff^-norm by em^ (U) thanks to (|2.3ip and Corollarv 12.201 
We deduce that Ra satisfies the estimate of the proposition. The fact that (|4.17p 
coincides with the identity given in the proposition is a consequence of the first 
point of Lemma 14.41 D 



We end this section with another linearization formula that describes the 
way 1/H_ g^^p,-. departs from g^/j(^a) (when a = 0, both terms are equal). 
This proposition will be used to express tpr^) in terms of the "good unknowns" 
C, and ■0(o,) in the proof of Proposition 14.91 In this sense, it can be viewed as 
a generalization of the transmission problem solved in Proposition 12.171 The 
additional term with respect to the case a = is responsible for the Kelvin- 
Helmholtz instabilities, as shown in the next section — of course, in the water 
waves case {p~ = 0), one has g''^ = g and ip'i^ \ — ipia) and this destabilizing 
term disappears. 



Proposition 4.5. Under the same assumptions as in Proposition \4.S\ one has 

^?^^U=^t^ (|a|<iV-l), 

-^c^±Vf„) = g^o^) - P^epiH^g o {g^r^^ ■ [d^av^D + r^, (l«l = N) 

where r^ satisfies the estimate, with |CI<iV+i/2> o,s in |^. jg[ ), 

Remark 4.6. A consequence of this proposition is that for all \a\ < iV, 

P*IWf„)|2 < m^(f/)(l + eM'/VVlCl<iv+i/2>). 

Proof. As for the proof of Proposition 14.31 we focus on that case |a| — N. 
Applying d" to the identity -^Q'^ij}'^ — g^ and using (|2.30p and Lemma H^ i 
we obtain 

where r^ is a sum of terms of the form, 

Bts,. = d^g{d'\: ■ ■ .:d^\)d'^ - -^d^g^{d^\, . . .,d^'C)d'^P^ 

with X;Li ki + |i5| = TV , < |(5| < TV - 1 and |i'| < TV. The formula of the 
proposition is obtained after replacing T [U] by its explicit expression (|4.ip . 
We now must control p^{g^)~^B^g ^ in i?^' (R''). The + case is more difficult 
than the — case because p~ is possibly very small while we always have p"*" > 1/2. 
Therefore, we only show how to control -B^^ j moreover, only three endpoint 
cases deserve a proof, namely: 
(i) The case j = |t^| = 1 and \6\ = N — 1. One then has 

B^s,. = dg{d^\)d'i,-^dg+{d^\)d'i;+, 

= ^dg+id^\)iJ-'d'^p - aV) - -^g+J-'dJ{d'\)j-'d'^, 

where J — J[Q and where we used the identity g — -^g^JlQ^^. Since by 
(|2:38l) . -0+ = J[C]"V, we thus get ~ 

B+s. = -^dg+{d^'OJ-'[J,d']J-'ij~^g+J-^dj{d^\)j-'d'ij 

-= B1 + B2. 
We now need the following lemma. 
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Lemma 4.7. For all 5 G N''+\ \5\ < N - I, one has 

|[J,5>|^3/2 < £p-m~(C/)(l + Ai^/4|C|<A.+i/2>). 

Proof of the lemma. Recalling that J[C] — p^ — Ptpf^Q)^^Q^i ^-nd writ- 
ing V = {g^Y^Q'^u, we get {d^ ,J\u = p^ ^{d^ ,{g^Y^ g^\ and we observe 
moreover that 

Since both terms in the r.h.s. can be handled similarly, we just give the details 
for the first one. It is a sum of terms of the form 

c,,8'.. = {g-r^d^g-{d'"c,,...,d^\)d''v, 

with X]i=i kN + 1^1 = 1*^1 ^^d \5'\ < \S\. We can control all these terms thanks 
to (|2.35p . (12.361) and Remark 12.161 except the most singular one (with respect 
to C), which corresponds to j ^ 1, l^ ~ S and 6' — 0. For this term, we write, 
using (ESDI), 

Using the inequality \d^(\^3/2 < A^~^^^|CI<a'+i/2> to control the first term and 
Remark l2.10l to control the second one, we get 

< £m^(C/)(l+//4|C|<^+l/2>), 

where the second inequality is a consequence of Corollarv l2.20l The lemma then 
follows easily. D 

Using (|2.16p and the lemma, we get the desired control on Bi ; for B2 there 
is no particular difficulty and we omit the proof. 

(ii) The case j = |^| = 1 and |t^| = A^ — 1. One proceeds exactly as in the proof 
of Lemma 14.71 
(Hi) The case \d\ ~ 0. j = 2 and |t^| = A^ — 1. One then has 

B+,^, = d'gid^\,d^\)i;-^d'g+{d^\,d^\Y+, 

and the scenario is the same as in case (iii) in the proof of Proposition l4.3l D 
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4.3 The quasilinear system 

We show in Proposition 14.91 below that the two-fluid equations (|1.19l) can be 
"quasilinearized" . In these quasihnearized equations, an operator plays a cen- 
tral role with respect to the formation of Rayleigh-Taylor/Kelvin-Helmholtz 
instabilities. We thus call it the instability operator. 



N 

(T,T 



Definition 4.8 (Instability operator). Let to > d/2, T > and (CV') S E^ 
{N > io + 2) be such that (|2.ip is uniformly satisfied on [0,r]. The instability 
operator Ins[U] is defined as 

Ins[U]f = af-e'^Ip+p-lV^j • £[C](/[y±]) - ^V • ^ [eV7^VC]V/, 
with 

For notational convenience, we also introduce the matricial operators 
and we also define Ca \U] as 



r \jn -( ° "^^(" 



(4.22) 



We can now state the main result of this section. 



Proposition 4.9. Let T > 0, tg > d/2 and N be as in Proposition \4.3\ If 
U — (^, ijj) E E^T satisfies 12. 1\) uniformly on [0, T] and solves fl.l9\) then for 
all a E N'^+^ with I < |a| < N, the couple [/(„) — (C(Q)iV'(a)) solves 

dtU^^-, + A [;7]f/(„) = e(i?„, S^f, (|a| < TV), 

where Uia) — {C(a) i i' {&)Y' '^''^d. the residuals Ra and Sa satisfy the estimates 

\R»\Hi + \S»\Hl^2<m^iU){l + ep'/'p+p-\lV^JUCUN+i/2>), 
uniformly on [0,r]. 

Remark 4.10. The "quasilinear" system we refer to is the system of evolution 
equations on (C/(Q))o<|Q|<Af (with C/(o) — U) formed by (|1.19p for a = and the 
equations of the proposition for 1 < |a| < A^. 

Proof. Throughout this proof, we denote by (f)^ the solution of the transmis- 
sion problem \2. 36]) provided by Proposition \2.17\ and write i/j^ = 4>^\^_„ (or 
equivalently, we define ip^ by 112.38]) ). 
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The first equation is obtained directly by applying 9" to the first equation 
of (jl.igp and using Proposition |431 

For the second equation, we focus on the most difficult case, namely, |a| — N. 
We decompose first a = (3 + j £ N'^^^, with I7I = 1. Inspired by [35], we apply 
first d'* to the second equation of (|1.19p and use the definition (|4.2p of V_ and 
wj^ to get 

- -lP^w^d\gm = -^^d^^{e^O- (4.23) 

We are now going to apply d^ to this equation. We will thus find many terms of 
the form d^ ifo)', in order to check that most of them can be put in the residual 
term (i.e. thail^"! |(9^(/g)| /,i/2 < em^(C/)), we rely heavily on the product 
estimate (|4.16p . For the sake of clarity, we also introduce the notation 

a^b ^^ \a-b\^^/.<em''{U){l + e^i^/^p+p-\lV^lU\(:\<N+l/2>)■ (4.24) 

Lemma 4.11. Under the assumptions of the proposition, and with a ^ j3 + ^ 
(\'y\ — \), the following identities hold: 



e 



^P^V"^ ■{d^^,ur^}\/d'-'C, 



-dl^ip^ur^ d^ig^P)) ~ -P^ur^d''{gij) + -p^{d^,ur^}d^{gi;), 
M M- P- 

where {d^ ,w^} = X]-i=i PjdjW^d^\ with J3^ + jSjGj = (3, is the usual Poisson 
bracket. 

Proof. Let us address the first assertion of the lemma. Since ep^d^(V_ ■ 
{yd'^ip^ — ewj^Vd'^C)^ is a sum of terms of the form ^^'.^" = ep^d^ V_ ■ 
dP" {Vd^i}^ - ew^Vd'''C), with (3' + /?" = /3, and remarking that 

^o,/3 ~ e^V^ ■ VVf„) + eV(Z'^ ■ Vw±)a"C - ^VZ"^ • {d^ ,w^}^d'^C. 

we are led to prove that Api ,pn ^ if |/3"| < A^ — 1. The most difficult con- 
figuration corresponds to /3" = or |/3"| = A^ — 2, and we thus omit the proof 
for the intermediate cases. Moreover, since the proof follows the same lines for 
(3" = and |/3"| = iV — 2, we just give it in the latter case. By (14. 16^ . we get 

|A/3',/3"|^y^ < £P^|a^V*lff'o+i(|a^".9|2 + \d^"9\Hlp), 

with g = Wd'^ijj^ — euf^Vd'^^. Since \f3'\ = 1, we deduce easily from Corollary 
12.201 that 19'^ l^^|_f/to+i < ^^(U); moreover, d^ ' g can be put under the form 



"The notation A < m" {U) is used as a shortcut for A{t) < m"{U{t)), for all < t < T 
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Since the last two components of the r.h.s. are bounded from above in LF' and 
H^' by m^ {U) (this is a consequence of (I4.16|) . Corollarv l2.20l and the identity 
l/ljji/2 < l/ljji), we have obtained that 

< em~(C/)(l+ Y. E^KUl'^) 

\k\<N 

< em^iU) (1 + efi'/'p+p- \ IZ*] U |CI<^+i/2>) , 

where Remark 14.61 has been used to derive the last inequality. 
In order to prove the second assertion of the lemma, we proceed as above to 
check that it is sufficient to prove that ^p±5'''w±9'^"+'^(^'0) -- if /3' + /3" = /3 
and |/3"| < N ~ 2. We give the proof of the most difficult case, corresponding 
here to /3' = /?, /?" = 0. Thanks to dUS]), Corollaries [2201 and [221] and Lemma 
14. 4[ it is enough to prove that 

^laVl^y. <m^(C/)(l + eMi/Vp-|lZ^lloo|Cl<iv+i/2>). (4.25) 

Recalling that wj^ = 'W^[(^]ip and differentiating this relation with respect to 
d" yields after multiplication by e (in order to use the convenient notation ^) 

(1 + eVlVCH— a'^w* + 2e2^V9^C • (ewz^VC) 



Since moreover, one can deduce from ()2.30p . p.3ip and (|2.33p that 



M+ " ' 7^+-' '^^^ 



,±^ 



we get the formula (after multiplication by p 

(1 + eV|VCP)^aV ^ ^(^±(p±^± p + £2^VC • v(p±^±p. 

Dividing this formula by e, it is then easy to deduce from Remark 14.61 (and 
(|2?20l) to handle the ffist term of the r.h.s.) that (1425)) holds. D 

Lemma 14.111 allows one to put most of the terms of the second equation of 
9'^ (|4.23p in the residual. However, this is not the case of the Poisson brackets 
that appear in both identities of Lemma 14.111 because they are not uniformly 
controlled by the energy, since they require (when evaluated in _ff^ -norm) a 
control of {N + 1/2) derivatives of C- Using the control of the derivatives of 
order (N + 1) provided by |CI/fi." would induce a singularity when the surface 
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tension goes to zero. Fortunately, this singularity disappears when one takes the 
difference of the expressions considered in Lemma 14.111 which is the quantity 
that appears in 5'^ (|4.23p . Indeed, one has 

p- 

^ 0, 

where (|2.30p has been used to derive the second identity. 

It follows from this analysis that (?'^ (|4.23l) can be written under the form 

/i - l3o e^/JL 



Using the first equation, one can replace —Qi! by dtC, in this expression, 

9t^(„) + a9"C + elp^V^ ■ V^± )1 ^ -i--l=9"^(£VMC), 
with o as in Definition 14.81 We can then use the identity 

[/Vl = (/^>1.9^1 + 1/^1(5^) 
to obtain 

9«V(a)+a9"C + g(Z^)-VW)+4^^I-V(p^V^f,)) - -^ ^°^J^^^ - (4-26) 

The quantity V(p^?/'r^x) depends on ijj/a) ^-nd d"C ^^id on their first order deriva- 
tives. This dependence is trivial in the water waves case since {p^ipia)) ~ 2^(")' 
when p~ y^ 0, it is responsible for the Kelvin-Helmholtz instabilities. The fol- 
lowing lemma makes this dependence precise. 



Lemma 4.12. One has, with £ [Q as defined in |^.P[ j, 

-e'pp+p-lV^j-^KMV^WO- 
Proof. From Proposition 14.51 we get 
1 



j^^U = ig^)-'gy^ia)-P^s^,K^ig^)-'g{g^)-'v■{^'^av^B+ig'- 



We deduce from this expression that 

elZ^l • V(p±Vf„)) = f IZ^I • V[(l + 2p-H-ig-)-' o (^)V'(„)] 
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The lemma follows from the fact that e[Z=^l • ^{p^H^{g^)~'^r^) ~ 0, which 
is itself a direct consequence of the bounds on r^ given in Proposition 14.51 D 

Thanks to the lemma, we can write (|4.26p under the form 

5,V(a) - el [C/]*V(„) + aC(„) - e'/ip+p-[Z±] • E [C](C(a)Ii:*l) 

^ ^ --d^l^ie^O- 



Boe.^ 
The proposition is thus a direct consequence of the observation that 



Bo Ey/JI Bo Bo 



D 



5 Main results 

5.1 The stability criterion 

5.1.1 A first criterion 

We state and comment in this section the stability criterion that ensures the 
existence of a "stable" solution, that is, of a solution that exists on a time scale 
consistent with physical observations. 

Before stating this criterion, let us recall that e(C) is defined in Proposition 
IMl as 

(g-^v-vy-v) 

with g = p^jY+g^ — p^-^g^ . We also define c(C) as 
c(C) = e(C)^(l+eV|VClL)^/^ 

and recaU that T = {p+ p-f^iP-±P-}^. 
- - if 2 _ 4(7 ^ 
The (nondimensionalized) stability criterion can then be stated as 

Tc(C)max|a«|y±l|'' <info, (5.1) 

with a = 1 + e|p* {pt + eVj^ ■ V)w*]] and where we recall that 9" can be either 
a time or space derivative. 

In order to measure by which extent this criterion is satisfied, we introduce the 
function c)(C/) defined as 

£)([/) =infa-Tc(C)max|a"|l^=^l|'* ; (5.2) 

quite intuitively, we expect the solutions to be more "stable" for large values of 
J)([/) than for smaller ones. 
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Remark 5.1. In the flat case C = 0; one can check that 



c(0) — sup 



X 



x>o (1 + x){p- tanh(iJ+j:) + p+ tanh(iJ x)) 

In the general case, it is possible to obtain an upper bound on c(C) by tracking 
the constants in Proposition l4.2l The most important information, however, is 
that c(C) ^ 1 for all the possible values of e and /i. 

5.1.2 An alternative version 



We show in Theorem 15.61 below that (|5.ip ensures the existence of a "stable" 
solution if the fluid depths do not vanish. If we assume moreover that the 
jump of the horizontal velocity is not identically zero at i = (which is the 
general configuration) , it is possible to give a more elegant version of the stability 
criterion, namely, 

Tc(C)|lZ^Hl<infa, eH) 

I loo gi 

Consequently, one replaces c)(C/) in (15.21) by 

y([/) = infa-Tc(c)|ii:^]|l. mi) 

gd I loo 

Remark 5.2. It is of course possible to give a version with dimensions of the 
stability criterion ()5.1[|) . Remarking that {p'^ + p^)g'a — [9^P* U=cl' where 
P"^ is the pressure in the Euler equations (11.41) . we obtain 



l^=«^ 4a(p^ 



-d.P^u-J>i-?^^r^^^AO 



w 



where u — [y* U=cl- This is the criterion stated in (jl.ip in the introduction. 

Remark 5.3. In the water waves case {p~ = 0), (|5.1|) and (15. lip reduce to the 
condition inf^d a > (or to the Taylor condition —dzP\^^, > in physical vari- 
ables). This condition is the same with or without surface tension^^l. This con- 
dition is usually not imposed in the literature when surface tension is presenio; 
the reason why it appears here is that we are only interested in "stable" solu- 
tions. In the water waves case, this amounts to discarding configurations where 
the water is above the air; though there exists a local solution thanks to surface 
tension, it is obviously not "stable" in the intuitive meaning of this word. 

5.1.3 A practical criterion 

We end this section with the derivation of a "practical criterion" that can be 
used in (shallow water) applications to have an a priori idea of the stability of 
an interfacial wave if we are given its physical characteristics (amplitude and 



^^with the convention that (p )^/(T = if p =0 and a = 

^■^but it is of course imposed in the references considering the zero surface tension limit, e.g. 

fHllliollil] 



56 



wavelength) . 

The nondimensionahzation performed in Appendix |^ (and supported by the 
experimental data of [37] for instance) is such that V_ is roughly of size one 
(this remains true as long as the dynamics does not depart too much from the 
linear analysis, which is the case for "stable" waves until their breaking point). 
It follows that a is roughly of size one and that a practical criterion can be 
stated as 

T <C 1 : Stable configuration; T ^ 1 : Unstable configuration. (5.3) 

When T '-^ 1, a refined analysis is of course needed and the full criterion (|5.ip 
should be used. 



Remark 5.4. It seems that the wavelength does not play any role in (|5.3p . Its 

umpti 
is only true in shallow water, i.e. when H^/X^ <C 1. 



contribution is however hidden in the assumption that V_ is of size one, which 



5.2 Initial conditions 

The main step in the proof of the existence of stable solutions is to prove that 
the energy £^{U) defined in (I4.10[) . namely, 

C(t^) = |VV|?,*o+^ + E IC(a)|?,i+|V'(o)l|l/., 

aeN''+MQ|<Af 

is controlled for positive times by the energy at t = 0. Since £^{U) involves 
time derivatives, we have to specify which sense we give to the initial energy 
£^iU°). If we denote, for all a € N-^+i, a = {ao,ai, . . . ,ad) so that d"" 
corresponds to the time derivatives of 9" , the problem is to choose initial values 
^{a) fo'^ iU{a))\t=o ('^ith U(a) = (C(a) , V'(q) )) whcu oq > 0, in terms of U° and 
its spacial derivatives. 

This can be done by a finite induction. When ao = (no time derivative), we 
take of course 

U^^^ = (a"C°, 9"^° - eu/^d^Cf with u;0 = [p±w±[£C"]V'°l- 

Now, let 1 < n < A^ and assume that C/(^) i^^^ — U9p^ has been chosen for all 
/3 £ N'^"'""'^ with /?o < n. We remark that for all a with ao ~ n we have 

^(") \t=o = {dtC(a'),dt^(a') + sdtwd"' C) ^^^^, 

with a' = (ao — 1, ai, . . . , ad), and we are therefore led to set initial conditions 
for 9fC/(Q/), which is achieved by using Proposition 14.91 

The initial energy, which we denote slightly abusively by £^ {U^) is therefore 
defined as 

£,^([/°) = |V^°|?,.„,. + Y. IC(%l?.i + l<a)l^v., (5.4) 

Q,gN<i+i,|a|<A' 

with U9 N as constructed above. 

(a) 
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Remark 5.5. One proceeds in the same way to give a formulation of the stabihty 
criterions (|5.ip or (|5.1[p at time i = in terms of C/° and its spacial derivatives. 

5.3 Local existence of "stable" solutions 

We state here a theorem ensuring that the interfacial waves equations (|1.19l) 

admit "stable" solutions if the depth of both fluid layers does not vanish for 

the initial condition, and if the stability criterion (|5.ip or (|5.1[p is satisfied. By 

"stable" solution, we mean two things. 

Firstly, the existence time must not shrink to as tr — >■ 0. More precisely, the 

existence time is not measured by the size of the surface tension coefficient a 

but by the size of J)(C/). 

Secondly, the existence time must not vanish for "acceptable" values of the 

physical parameters e and /i, namely < /x < 1 and < e < 1. 

Recalling that the quantity 8^ {U'^) has been defined in (|5.4p and E^rp in 
(|4.13p . we can now state the following theorem (whose proof is given in ii5.5p . 

Theorem 5.6. Let to = 3/2, N > 5 and U° = (C°,'0°)^ & L^iR'^) x H^^^{R'^) 
be such that 8^ [U^) < oo and that satisfies the nonvanishing depth condition 



3/i^„,>0, ^inf^(l±£±C'(X))>/. 



// moreover U^ satisfies the stability criterion i5.1\) then there exists T > and 
a unique solution U G E^rp to il.l9]) with initial condition U^ . Moreover, 

Vie[0,T], 8^{U{t))<CiT) 



re- 



and T depends only on Tn^([/°) and '^{U ) (as defined in {4-14^ o,i^d \5.^) 
spectively). 

If the jump of the horizontal component of the initial velocity is nonzero, 
one can replace the stability criterion (|5.1I) by its alternative version (jS.Hp . and 
get the following theorem. 

Theoremimr. Let to = 3/2, N > 5 and U" ^ (C°,0^ e L^{R'^) x H^^^{R'^) 
be such that 8^ (U^) < oo and that satisfies the nonvanishing depth condition 



3/ii^„>0, mf^^il±e^C\X))>h: 



together with the condition cu^^ = [Hi 7^ 0. // moreover U^ satisfies the 
stability criterion h5.1\]) then the conclusion of Theorem 1 5. 61 still holds, with T 

depending on m iU^), iU ) and -r-^, • 

Remark 5.7. 

(1) The existence time also depends implicitly ^^; since we address here physical 
configurations where H^ and H^ are of the same ordea^i. this dependence is 
harmless. 



^^As said in the introduction, this assumption is made for the sake of clarity and the method 
presented here could be adapted to other configurations 
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(2) There is also an implicit dependence of T on 1/Bo. Since we are not interested 
in situations where capillary forces dominate the effect of gravity (e.g. droplets), 
the Bond number is always larger than one and the dependence on 1/Bo is 
harmless. 

(3) It is straightforward to deduce from Theorem l5.6l some convergence results as 
p~ and/or a go to zero. For instance, convergence to the water waves equations 
without surface tension is obtained when p~ and a go to zero with T ^- - and 
in particular if {p~)'^ /a — > 0, a result established in 46 for a liquid droplet with 
rotational effects under the slightly stronger condition (p^)^ < cr^/^. Theorem 

l5.6l shows moreover that it is possible to take the shallow water limit at the same 
time, which justifies the use of the standard (one fluid) shallow water models 
for the air/water interface (see Section 16.11 for more details) . 

5.4 Persistence of "stable" solutions over large times 

We have seen in Remark [5.3l that in the water waves case without surface tension, 
the stability criterion (|5.ip reduces to the Ray leigh- Taylor criterion inf^d o > 0. 
If this condition is satisfied, we know from [3] that the solution U exists for a 
time scale 0(l/e). When e ^ 1, this results coincides with Theorem 15.61 but 
when e <?C 1, Theorem 15. 6l is weaker than [3j since it provides an existence time 
of order 0(1) only. 

We propose in this section a stronger version of the stability criterion (|5.ip 
(one can adapt (|5.1[|) in the same way) that allows us to show that the solution 
furnished by Theorem 15.61 persists over larger times when e <^ 1. This stronger 
criterion can be stated as 

e-^''T^c(C)max\d°'lV^jf <mia, with 0<7<1 (5.5) 



(in the case 7 = 0, (|5.5|) coincides of course with ()5.1|) ). 

Theorem 5.8. Let to = 3/2, N > 5 and U° ^ iC°,i^°V G L^i^'^) x H^/^ 
be such that £^ [U^) < oo and that satisfies the nonvanishing depth condition 

^hi,„>0, inf (l±£±C°W)>/ii^„- 

// moreover U'^ satisfies the strong stability criterion \5. 5\) for some 7 G [0, 1], 
then there exists T > and a unique solution U G E^ ^-^rp to lil.19]) with initial 
condition U^ . Moreover, 

yte[0,e-''T], £^{U{t))<C{T) 

and T depends only on m^(C/*') and d{U ). 

Remark 5.9. If (|5.5p is satisfied with 7 = 1, we recover the same time scale 
as the one provided by [3 for water waves. We refer to Section |6] for further 
discussion. 

Remark 5.10. The practical criterion corresponding to (|5.5p is obtained by re- 
placing T by e~'^^T in (j5.3|) . 
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5.5 Proof of Theorems ^M ESf and ISTSl 

We give here the proof of Theorems 15.61 and 15.81 The proof of Theorem | 
requires only minor adaptations with respect to the proof of Theorem 15.61 we 
briefly indicate them in footnotes. 

5.5.1 Preliminary results 

We construct below a symmetrizer for the quasilinear system (|5.8I) . It is of cru- 
cial importance that the energy norm associated to this symmetrizer (see (|5.9p 
below) be (uniformly) equivalent to the energy £^ {U) previously introduced in 
(|4.10l) . A first ingredient is the coercivity property oi -Q established in Propo- 
sition [2211 The second ingredient is the following lemmeothat shows that, up 
to lower order terms that can be controlled and under the stability condition 
(|5.ip . {lns\U]u,u) is uniformly equivalent to |u|^i. 



Lemma 5.11. Let T > 0, to > d/2 and U = (C, V') € E^j, be such that ![2l\) 
and i5.1]) are uniformly satisfied on [0, T] . Then the following identities hold 
uniformly on [0,T], 

(1) For all u G H^{W^), and with xn}{U) as defined in ^.14\ ), one has 

{lns[U]u,u) <m\U)\u\jji. 

(2) There exists a numerical constant C2 such that 

\KU)\u\hI < {lns[U]u,u) + aiU)\u\l+b{U)\u\l^,,,, 



|2 



wttha{U)^e^p+p-c{0\l}ri\^ andb{U)=e^^p+p-C2c{0\l}ri\H.„+i, and 
d{U) as in (El). 

(3) With I • 1(1/2) as defined in \4-.12^ , one has 

e'^p+p- max |9"C*l|lk|^i/2) < m^UMni. 

Proof. By definition of Ins[U] (see Definition 14. 8p . one has 

{lns[U]u,u) = {au,u)~e'pp+p-i^{ulV^l),ulV^i) + J-(:?c Vu, Vu). 

- - Bo 

In order to give some control on the three components of this expression, we 
need the following straightforward estimates, 



(infa)|u|^ < {au,u) < m\U)\ 



U\: 



Vu\l < (5C[ev^C]Vu,Vu) < \Vu\l. 



(1 + £V|VCIL)^/' 



^''if one replaces l|5.1|l by 1)5. lip then the estimate of the third point only holds for a = 0. 
In the case |«| = 1, we must therefore replace the estimate of the lemma by 
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and, using ProDOsition l4.21 



< {^{ulV^}),ulV^J) < 'Mui + ^\D\y/^{u 



^-^J)|^ 



where we used the estimate 

V/'lE^lloo||i5|i/2«|^ (5.6) 

to derive the last inequahty (here C2 is a numerical constant). 

The upper bound on {lns[U]u, u) given in the lemma is a direct consequence of 

these inequalities. For the lower bound, remark that they also yield 

{lns[U]u,u) >A- a{U)\u\l - 6(f/)|w|^_i/,, 

with 

1 iVwP 

A = (i„f a)|„|^ - sfT^(v)Hln + g^ ^^ _^ evivk)'^' 



Basic calculus shows that under the stability condition (|5.ip . one has 2 A > 
3(C/)|u|^i, and the result follows. Remarking that A remains positive if one 

replaces \IV_ ]|oo by max|Q,|<i |9"|il'''|| in a{U), one gets the last point of 
lemma as the first one. D 

5.5.2 The mollified quasilinear system 

Let X : H^ ~^ US. be a smooth, compactly supported function equal to one in 
a neighborhood of the origin. For all i > 0, we denote by J^ the moUifier 
J'^ = x(''l-C'|). Let us consider the following mollified version of the two-fluid 
wave equations (|1.19p . 

M 

dti^ + JX + e\r{lp^\W^^\'j-^{l+e'n\^(:mp^{w^n) (5.7) 



Bo e^/JI 



From standard results on ODEs, (|5.7I) has a unique maximal solution [/' = 
(0''') with initial data (C°, "0°) on a time interval [0, T^^^). Proceeding exactly 
as for the proof of Proposition 14.91 one can check that for all a G N'^"'"^, 1 < 
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|a| < TV, the quantity C/(^„) = (((„)>(„)) (with Q„) = d'^C and V^„) = 9«i^' - 
eIp±w±[eC']V'i5"C) solveS 

att/(„)+J^^[(7]C/(„) + J'S[C/][/(„)+J^C«[(7]C/(e,) =e{.rR^,rSo, + S'^f (5.8) 

(as in Proposition I4.9[ the operators J? [[/] and Ca[t^] niust be removed if \a\ < 
N), where Ra and Sa satisfy the same estimates as in Proposition 14.91 and 
where S'^ is given by 



S'^ = -(l-r)[lp^dtw^ 



*)J 



Our strategy is to derive energy estimates for the system of evolution equations 
on {U^.)o<\a\<N formed by (|5.7p and (|5.8p . After showing that these energy 
estimates do not depend on the molhfying parameter l, we wiU construct by 
compactness a solution {U(^a))o<\a\<N to the limit system corresponding to t = 
(that is, J' = 1) on a time interval [0,rCT]. The component corresponding to 
a = of this solution furnishes a solution to (|1.19l) . Without further analysis, 
one has a priori T^ ^> as cr — > 0. An energy estimate on the exact equations - 
as opposed to the mollified ones ~ allows us to prove that the existence time can 
be taken as in the statement of Theorem 15.61 We then prove that the solution 
persists over larger time if the stronger criterion (j5.5|) is satisfied. 

5.5.3 Symmetrizer and energy 

Since d — 1,2, we can take tg = 3/2 > d/2 and any N > 5 satisfies the condition 
Ji4.15}^ . A symmetrizer for the "quasilinear" system (|5.8p ,v is given by 

S[U] ^ S^[U] + s^iu] + s'^iu], 

where 5-'[C/] (1 < J < 3) are diagonal operators defined as 

S'[U] = dil,gilns[U],-g), 

S^[U] = diag(^3C(a)[£\/7^VC],-^(„)), {orS^[U]=0 if \a\ < N), 

S^[U] ^ diag(a(C/)+6(C/)A-\0), 

where the constants a{U) and b{U) are as in Lemma [5. Ill 

Multiplying (|5.8p by 5 ^ [U] symmetrizes the highest order terms, but since S ^ [U] 

is a second order operator, the subprincipal terms need to be symmetrized 

also (otherwise the commutator term is not controlled by the energy) and this 

explains the presence of 5 ^ [U] . The S ^ component contains two terms of order 

and —1 respectively; its role is to ensure that {s[U]U,U) controls the energy 

C(C/). 



-•^^for the sake of clarity, wc write 0^\ and ipict) instead of f .' . and ip'j . when no confusion 
is possible. 
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This discussion motivates the introduction of the energies Tl{U), for all 
1<1 <N, 

tUu)^ V T:iU), (5.9) 



o<|q|<; 



with 



2 1 (5.10) 

J-^(C/) - |C|2,, +-(V,,(^^[0]^) (a = 0) 

M 

where [/(„) = (C(q)j V'(q))i ^nd This definition raises at least two important 
remarks: 

1. The component 5^[t/] does not appear in the energy; the reason is that 
(5^[t/]9tC/(Q),f7(a)) can be controlled by the energy without using the 
equation to replace the time derivative by space derivatives. This is made 
possible by the fact that the energy controls both space and time deriva- 
tives. This method was introduced to handle surface tension in the water 
waves case in [52] (see also [10]) and can be adapted here with small mod- 
ifications. 

2. The presence oi S^{U) induces nonsymmetric terms. Since these terms 
are of lower order, they do not raise any difRculty as far as local existence 
is concerned. However, they may destroy the stability criterion since they 
may interfere with the zero order component (namely, a) of Ins[U] - or 
at least restrict its range of validity (very small values of p^ for instance) . 
This explains why the zero order term of S^{U) is treated with special 
care in the computations below. 

Let us finally remark that, owing to Proposition l2 . 24l Lemma [4.4l and Lemma 
15.111 the energy F^ {U) is equivalent to £^{U), as defined in (I4.10p . 

C(f/) < (M + ^)^^iU) and J-f ([/) < m\U)£^iU), (5.11) 



with i){U) as in (lO)) . 

5.5.4 Energy estimates 

For notational simplicity, we write E = £[C], ^ = !!(,[£ y^C], T — 1 [U] etc. 
We first consider the case a ^ 0. Taking the L^-scalar product of (|5.8p with 
(5i+53)[/(„)+52C/(a) yields 

((5^ -l-53)5t [/(„),[/(„)) -1- (5,f/(„),52[/^^^) 
= e(J'(i?c.,5„ + 5;.)^,(5i+53)C/(„)+5^C/(a)), 
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where we used the fact that (j''JiU(^a),S^U(^a)) — and (J'^-^C/(q),5qC/(q)) + 
{j''CaU^a),S^U(^a)) ~ 0. By definition oi F^{U) we get therefore 

, 7 

-(j-,"(C/) + (C/(„),5^{/^^^)) ^^A,+i3i+i?2, (5.12) 



with 






As = -(j'CaC/(a>,5^C/(a)), Ae = -((5^ +5^)J'®t/(„),C/(„)), 

and 

We now give some control on the difi^erent components of the r.h.s. of (|5.12p . 

- Control of Ai. We can decompose Ai = An + A12 + ^13 + A14 with 

All = {[dt,Ins[U]]C{a),C{a)), M2 = ij[dt,g]^p{a),^(a)), 

One gets easily that \Aij\ < em^{U) for 2 < j < 4 (this is straightforward for 
A13 and Ai4, and a consequence of Proposition 3.6 of [3] for A12). We thus 
focus our attention on An. First remark that 

An = ((atQ)C(a),C(a)) + ^([at,i*C]VC(„),VC(a)) 

~ 2eVpV(2:(C(a)li:*l),C(a)at[z*i) 

It is straightforward that the first two components are bounded from above by 
em^([/). To control the other two components of this identity, we use Proposi- 
tion 22] and (|5.6p to obtain 

\An\ < em''{U){l + e'^p+p-\lV^j\l\\D\^%^^\l) 

+e^VJlp+p-\lV^\oomV^\oo\\D\'/\^^)\l. (5.13) 

Thanks to the third point of Lemma [5. Ill we get thereforq^ I An] < m^ (U). 

- Control of A^. From the definition of 5^ and the fact that time derivatives 



^^if we assume US. lit rather than US. Il l then one has to adapt this estimate as indicated 
the footnote lfl4) l 
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are allowed in the energy (|4.10p . it is easy to get that |^2| < em^{U). 

- Control of A3. From the definition of 5^, it follows that 

A3 =a(C/)( J'- ^^(„),C(a)) +&(C^)(-^V'(a), J' A-'C(a)). 

Owing to p.23p . we get 

IA3I < MM|W(a)|2(a(C/)|*PC(a)|2+fe(C/)|*PC(a)|H-0 

and therefore, from the definition of a{U) and b{U), {A^l < em'^{U). 

- Control of A4 . Since one has 

and because |^(a)V''(a>|2 < em^([/) (see (|2.32p and Lemma W^ . we easily get 
that \Ai\ < em^{U). 

- Control of A5. As for A4, one gets without any difficulty that l^sj < em^{U). 

- Control of Aq. We decompose A^ = A^i + A^2 with 

^61 =- (5 V' S [/(„), t/(„)) and ^62 =- (5 V' S [/(„), L/(„)). 
We can first remark that 

ylei =-£(/n5[f/]J^TC(„),C(c))+e(-^^'1'*V'(a),V'(a))- 

We therefore deduce that \A^\\ < era^ (U) from the third point of Lemma [5. 11 1 
and the following identities, 

|(-^ J'T>(„),V(„))| < m^(L/),(5.15) 
M 

whose proofs are given in Appendix [C] The component Aq2 is much easier to 
handle and one gets readily that 1^621 < em^{U). 

- Control of Aj. One can write 

A7 = £(TC(a), ^!^(a)C(a)) " e(^>(a), -^(a)V'(d))- 

Using the first point of Proposition HTT] and the definition of ^(q) to control the 
first term, and the last point of Proposition 14. II and (|2.3ip - (|2.32p to control the 
second one, we get that \A-j\ < em^(C/). 

- Control of Bi . Thanks to the bounds established on Ra and Sa in Proposition 
and to the third point of Lemma lS.lli there is no pain to get |-Bi| < em^ {U). 

Control of B2 From the definition of 5^, (P?^ and (P3^ . we get that 



IS2I < m^(;7)(l + |C(a)kO < (1 + VBo)m^([/). 
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Gathering all the informations coming from the above estimates, we deduce 
from ((5l^ that for all 1 < |a| < A^, 

|(j-,"(t/') + ([/[„), 5^f7(^^)) < (l + ^/B^)m^(C/^). (5.16) 

For a = 0, we first rewrite the equations (|5.7p under the form 






(5.17) 



with 



AAi(t/) 



e/i 



(^mKIV'-^m 



A/'2(C/) = ^(lP^|VV^fl-i(l + eV|VCP)b±(m±)1 



eBo ymwvcF 



)vc. 



We then take the L^ product of ([5l7| with ((l - ^A)C, -^(^^[0]^) to get 



-J-°(t/) < e|AAi(C/)U.|CUi+e|^AA2(t/)|2|W|2 



< em''{U) 



(5.18) 



where the control of A/i(C/) and JV2{U) does not raise any particular difficulty. 
Summing (f5TT6| and ([gJ8| over all a e N'*+\ first with |a| < TV - 1 and 
then with |a| < N (and recalling that 5^ = if |a| < N), we deduce that 



-(J-;-^(C/')j < (l + VBo)m- ([/'), 
|(j-,^(C/^)+ 5] (L/(^„),5^C/(^))) < (l + VB^)m^(f/'), 

|Q|=Af 

and therefore, for any constant £, 

^?aiU') < (1 + \/B^)m^(C/'). (5.19) 

at 

with J^,^(f/^) = F^{W) + CJ-^^-Hf/^ + E|a| = iv(C^(a)''5aC^(^^.))- 

Now, it follows from the definition of 5^ that there exists a constant c^{U'') < M 
such that, 

\a\=N 
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As long as £ > C3{U') (and C < M), one has therefore, 

With the help of (|5.1ip . this in turn implies that 

9 ~ 



C(^') < 2(A/+ ^^)J-^(t/'); (5.20) 



it follows therefore from (|5.19l) that 
d 



^^-FJ^(C/')<(l + VB^)ft^(C/'), (5.21) 

where E^([/^) = C{M,^^m, i-, ^). 

5.5.5 Construction of a solution for times of order 0(1) 

We deduce classically from (|5.2ip that there exists T > as in the statement of 
the theorem such that j'^(C/'(i)) < C(T,T^{U°)) ior all t € [0, r*(l+VBo)"^] 
(we proceed as in >j5.2l to define the initial energy iP^ {U^)). 
Since this time interval does not depend on the mollifying parameter t, we 
can use a standard compactness argument (see e.g. [Tj for details in a related 
context) to construct an exact solution U G E^t ^^ (I1.19P - note that we use 
(|5.20p to deduce a bound on £^ {U) from the bound on F^ [U). 

As said in the introduction, the Bond number Bo is generally very large 
for the applications and the existence time of the solution constructed in the 
previous solution is much shorter than what can be physically observed. This is 
the reason why we paid a lot of attention to obtain energy estimates independent 
of Bo. In fact the presence of the prefactor (1 + vBo) in ()5.2ip only comes from 
the control of B2 in §5.5.41 Since B2 disappears if the energy estimates of ^ 35.5.41 
are performed on the exact equations (11.19^ rather than the mollified ones, one 
can replace (|5.2ip for the solution U constructed above by 

j^F^{U)<m^{U), (5.22) 

from which we deduce that it is possible to extend the time interval of existence 
to [0,r]. We omit the proof of the uniqueness, which is completely standard. 

5.5.6 Proof of Theorem [STSl 

In §5.5.41 all the Aj (1 < j < 7) and Bk {1 < k < 2) are bounded from above 
by em^(C/), except Ai and i?2. Since B2 is a mollifying error that vanishes for 
the exact solution, we can focus our attention on Ai, and more specifically on 
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All since An (2 < Z < 4) are bounded from above by em^{U). For An, we 
recall that we only have |Aii| < m^{U). More precisely, (|5.13p yields 

Using the third point of Lemma 15.111 the second component of the r.h.s. can 
be bounded from above by m^{U) but not by em^{U). The stronger stability 
criterion (j5.5p allows us to improve this estimate, as shown in the lemma below 
(whose proof is very similar to the proof of the third point of Lemma 15.111 and 
therefore omitted). 



Lemma 5.12. Let T > 0, to > d/2 and U = (C, V') € E^j. be such that l[2l\) 
and h5. 5\) are uniformly satisfied on [0, T] . Then one has 

e^^p+p-c{0 max \d'^iV^i\l\u\l/,^ < e^mi([/)|u|^i . 

Owing to this lemma, we get |Aii| < £"'m'^{U) and (|5.22p can therefore be 
improved into 

from which the theorem follows easily. 



6 Applications 

Physical implications of the stability criterion (|5.ip will be investigated in sepa- 
rate works, but in order to enforce its relevance, we describe briefly here some of 
them. We also show how Theorem l5.6l can be used to give a rigorous justification 
to two-fluid asymptotic models. 

There are mainly two different kinds of applications. The first one is pro- 
vided by air-water interfaces. In this case stability of interfacial waves is a 
consequence of the smallness of p^ (as in the zero density/zero surface tension 
limit considered in [46] for instance). The second application concerns internal 
waves for which both densities are comparable and stability is a consequence of 
the smallness of the shear velocity at the interface. 

6.1 Air- water interface in coastal oceanography 

We are interested here in applying the above theoretical results to typical con- 
figurations in coastal oceanography (see fll in the introduction). We want to 
investigate two things in particular: 



1. Can we use Theorem 15.81 to confirm that the standard models used in 
coastal oceanography (e.g. Boussinesq) do not need any correction due to 
the presence of the air? 

2. Can we still neglect the density of the air near the breaking point? 
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Before addressing these questions, let us give some numerical values for 
the physical parameters involved. The density of the air (at one atmosphere 
and around 20 Celsius) is p^ = 1.2kg • m^'^. For sea water, the density is 
p+ = 1025 kg- m"^. These values lead to p- = 0.001 and p+ = 0.999. For the 
air/water interface, the surface tension coefficient is cr == 0.073 N • m~^, and we 
take for the acceleration of gravity g = 9.81m • s^^. 

We also assume here that the air layer has the same height H^ = H^ = H 
as the water layer (this is of course not realistic, but as said in the introduction, 
the techniques used here could certainly be generalized to, say, an air layer of 
infinite depth, and the qualitative behavior should remain the samel ^1). 

6.1.1 Validity of the one-fluid long wave models 

We are interested in long waves models for which e+ ^ /i+ ^1. It is known 
[m [3] in the one-fiuid case that these models describe correctly the solutions of 
the water waves equations for times of order 0(l/e). 

Let us consider first a typical "long wave" over a depth _ff + = iJ = 5 m, of 
wavelength A = 35 m, and amplitude 0.1m. As said in Remark 15.91 in order to 
show that taking into account the density of the air does not yield significant 
modifications to the standard one-fluid theory, we have to check that the strong 
stability criterion (15.51) is satisfiecO, with 7 = 1. We rather check that its 
practical version (see Remark I5.10p is satisfied, which is the case since one 
computes easily that e~^T ^ 4.10^^^ <C 1. The solutions of the one-fluid and 
two-fluid models exist therefore on the same time scale. 

We also know from ^3 (one-fluid) and [11] (two-fluid) that these solutions 
are correctly described by standard long wave models (Boussinesq, KdV, etc.). 
The influence of the air density and the surface tension on these models being 
negligible (for instance, surface tension induces a modification of the dispersive 
term of the KdV model of less than 10"*%), this confirms the relevance of the 
one-fiuid models (without surface tension) to describe such phenemonons. 

6.1.2 Kelvin-Helmholtz instabilities for breaking waves and white 
caps 

A commonly observed value for e'^ near the breaking point of the wave is e"*" '-^ 
0.4, corresponding to an amplitude a = 6m for a depth H^ = 15m. One 
computes T r^ 0.27, which is not small enough to give any definitive conclusion 
on the presence or not of Kelvin-Helmholtz instabilities. We therefore have to 
use the exact criterion (|5.ip . This requires further experimental work, but it 
is likely that for some configurations, Kelvin-Helmholtz instabilities are at the 
origin of wave breaking. This might be the case for instance for spilling breakers 



^^Since p is very small, the influence of H on the typical length H is negligible. In our 
examples, only the third decimal of H is changed if H^ is taken equal to +00 

^*For these values, the Bond number is Bo ^ 1.7.10*. A local woU-posedness result relying 
only on the extra control furnished by surface tension, see i]5.5.5l would therefore give an 
existence about \/Bo£~^ ~ 6.10^ times smaller than Theorem 15.81 
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(as opposed to plunging breakers). It is indeed not understood which kind of 
singularity of the water waves equations creates such breakers; this might be 
because it is a two-fluid singularity similar to the "white caps" observed in 
presence of wind. The occurrence of these white caps is commonly explained by 
the linear Kelvin criterion stated in the introduction. It would be interesting to 
check whether a nonlinear version of our nonlinear stability criterion (15. ip with 
nonzero background explains part of the discrepancies of this theory (see |12|). 

6.2 Internal waves 

Waves at the interface of two immiscible fluids of comparable density are called 
internal waves. Waves at the interface of two layers of water of different density 
also carry the same name and are often described using the two-fluid formalism 
(see the recent review [IH])- The stability criterion (|5.ip is discussed here in 
these two configurations. 

6.2.1 The Koop-Butler experiment 

We consider here a classical experiment by Koop and Butler [35] . The upper 
fluid is deionized water {p~ — 998 kg • m~'^) and the lower fluid is Freon TF 
(p"*" = 1563 kg • m~'^). The depth of both layers in this experiment are H~^ = 
1.366 cm and H'^ = 6.948 cm. One gets therefore p~^ = 0.610, p~ = 0.390 and 
iJ = 1.989 cm. 

The interfacial tension coefficient a is not provided in |35| but the chemistry 
literature suggests that a = 0.005 N-m^^ is a reasonable value. Koop and Butler 
observed different interfacial waves of amplitude ranging from a = 0.034 cm to 
a — 0.68 cm. The dimensionless parameter T satisfies therefore 5.39 x 10~^ < 
T < 0.086. According to the practical criterion (|5.3I) , this correspond to stable 
configurations. Though the practical version of the strong stability criterion 
(|5.5p (see Remark l5.10p is satisfied, we do not discuss here large time existence 
because the viscosity of Freon TF should then be taken into account. 

6.2.2 The Grue et al. experiment and oceanic internal waves 

The configuration studied in [2 7) is slightly different from the experiment men- 
tioned in the previous section. Indeed, the authors do not consider two immisci- 
ble fluids but brine p"*" = 1022 kg-m^^ and water p^ — 999kg-m~^ (p+ = 0.506, 
p^ = 0.494). This configuration is thus closer to what is observed for oceanic 
internal waves. 

The difference with the case of two immiscible fluids is that there is a thin zone 
(pycnocline) in which the density is allowed to vary continuously between /o+ 
and p~ . This continuous stratiflcation is known to have a stabilizing role, which 
is important for small amplitude interfacial waves. 

It seems therefore that the physical framework studied in the present article 
differs from [27] . However, water-brine interfaces have been described with 
great success with models based the two- fluid equations (|1.3p - ([1.8p (see ^ and 
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references therein). It follows from our analysis that surface tension must be 
included to these equations if we want to use them to describe water-brine in- 
terfaces. This surface tension is somehow artificial and a natural question is to 
know which value it should be given. 

We propose here to use our analysis and the experiments of [37] to sketch a 
method to find a realistic value for a. The main idea is that its value should he 
such that the Kelvin- Helmholtz instabilities predicted by the theoretical analysis 
should coincide with those observed in the experiments. 

In [57], the depth of the two layers are respectively iJ+ = 0.62 m and 
H- = 0.15 m, so that H = 0.243 m. 

The authors report a set of measurements (Fig. 7 in [27] ) for intcrfacial waves 
of amplitude ranging from a — 0.033 m to a = 0.226 m. Kelvin-Hclmholtz insta- 
bilities are not observed up to a = 0.184 m but are present for the experiment 
with largest amplitude a — 0.226 m. It is therefore reasonable to assume that 
the critical amplitude is a = 0.2m. According to the practical criterion (|5.3p . 
this critical value should correspond to T ~ 1. We can then use this relation 
to get a = 0.095 N • m"-'^; this value is comparable to the air- water surface ten- 
sion, but slightly higher, which is not surprising since the stabilizing effects of 
the continuous stratification are taken into account with such a choice. Further 
investigation is required to check whether this approach can prove useful in the 
study of oceanic internal waves. 

6.3 Full justification of two-fluid asymptotic models 



We can use Theorem [5T8] to provide a rigorous justification of asymptotic models 
for internal waves, on the relevant time scale (to our knowledge, the only rigorous 
result so far is a justification of the Benjamin-Ono equations on the too-short 
time scale 0(Bo~^") that can be found in [43]): 

Many asymptotic regimes and models exist in the literature ([15] [20] [42] and 
references therein). Theorem l5.8l can be used to justify rigorously most of them 
along the procedure used in [3J for water waves models. 

• (Consistency) One proves that any family of smooth enough, uniformly 
bounded, solutions to p.l9p existing on the relevant time scale solves the 
asymptotic model up to a small residual. This is done for a wide class of 
regimes in [llj . 

• (Convergence) One proves that such exact solutions to the two-fluid equa- 
tions ()1.19|) remain close to exact solutions of the asymptotic model. From 
the previous step, this only requires energy estimates on the asymptotic 
model (generally much easier than for the full two- fluid equations). 

• (Existence) One proves that smooth, uniformly bounded family of solu- 
tions to (I1.19|) whose existence is assumed in the previous steps exist in- 
deed. As in the water waves case, this is the most difficult step. It is here 
ensured by Theorem 15.81 if the stability criterion (|5.5p is satisfled (which 
is the case for the stable waves of the above experimental observations). 
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Let us for instance implement this procedure in the "shahow-water/shalfow- 
water" regime characterized by e ~ 1, /u <C 1. In this regime, the interface 
elevation ( and v = Vtp are commonly described [13 [20] by the solution ((^°, v'^) 
of the following system (for notational simplicity, we give it for d = 1, but there 
is no additional difficulty to take into account the nonlocal generalization of this 
system to the case d = 2 derived in [TT] and studied in [25]). 



dtC + a. 






= 



dtv^+d,C+e-d, 



1 rp+h'{C?-p-h+{Cy 



a\2 



[p+h-{C'') + p-h+ic-w 



K) 



(6.1) 



with h-{C) = H^il - e-C) and h+{C) = H+{1 + e+C). We justify the 
approximation furnished by the solutions of this system in the following sense. 

Theorem 6.1. Let [/" = {C^,ip^)'^ be as in the assumptions of Theorem ] 5. 61 
and assume that 115. 1\) is satisfied uniformly with respect to p £ (0, 1). Let also 



V 



d^ip^ and assume that 



inf fl 



2 + - 

e p^p 



[H+^H-f 



[p+h-{C°) + p-h+{C')\ 



:{v'f) > (6.2) 



Then there exists T > such that for all p G (0, 1); 

(i) there exists a unique solution (C,t)") € C{[0,T]; H'^^^''^{ 

initial condition (C°7W°); 

(a) the solution U provided by Theorem 1 5. 61 exists on [0,T]. 

Moreover, with v — dxijj, one has 



\iC,v)-iC,vn\i 



([0,T] 



<R^) < pC{U°). 



'-) to {Op with 



Remark 6.2. The condition (|6.2p ensures the hyperbolicity of (16.1^ : it is not 
present for the standard one-fluid shallow-water equations. An asymptotic ex- 
pansion of the instability operator Ins [U] shows that the term responsible for the 
Kelvin-Helmholtz instability (the one involving E [C]) contributes to (|6.2p . which 
was suspected but not rigorously established so far. The shallow-water/shallow- 
water model (|6.ip is however certainly too unstable and the analysis performed 
here on the full Euler equations shows that the surface tension term should be 
kept to control high frequency instabilities. 

Proof. For the first point, we remark that |C°|^iv-i/2 -I- |i;°|^jv-i/2 < m^{U); 
the existence of a solution {Cv") S C([0,T]; 7?^-i/2(R)2) to dH]) with initial 
condition (C°,v°) (with T — T{m^{U))) is therefore straightforward because 
the assumption made in the statement of the theorem insures the hyperbolicity 
of (|6.ip . as shown in Theorem 1 of |28) . 

Taking a smaller T if necessary, it is possible to assume that the solution of 
Theorem 15.61 exists also on [0, T]; this theorem ensures also that supjQ j^j £a i^) 
is uniformly bounded with respect to p e (0,1). Thanks to the consistency 
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result provided by Theorem 4 of [TT], this ensures that {(,v) solves (|6.ip up to 
a residual /i(r^, r!^), with (r^, rj ) uniformly bounded with respect to /i G (0,1) 
in C([0,T];i?^(]R)^). A standard error estimate on the hyperbolic system (|6.ip 
yields therefore 

\{C,v) - {C,v'')\L^iio,T];miRr < mC(c/°), 
and the estimate given in the theorem follows from the continuous embedding 
7?i(M) cL°°(M). D 



A Nondimensionalization of the equations 

It is quite straightforward to nondimensionalize the space variable X and the 
interface deformation C using the quantities a, A, H^ and H introduced in !jl.3l 
We thus define 

A a 

where the tildes are used to denote dimensionless variables or unknowns. The 
nondimensionalization of the time variable t and the velocity potential tp requires 
the knowledge of a reference velocity c. Such a value can be obtained by direct 
observation, but is of course preferable to have a formula for c in terms of the 
above quantities. The best way to have access to such an information is through 
the linear analysis of the equations. Neglecting all the nonlinear terms in ()1.16p . 
one easily gets the system (see Remark 12.231 for the expression of Q [0] ) 

d,C-gm^ = 0. ^ith^[o] = |i^| t-h(F+|D|)tanh(i7-|i?|) 



5t^ + .g'C = 0, ^^J ' V+tanh(i7-|i:>|)+p-tanh(iJ+|L'|)' 

In the shallow water limit, that is, when H^ is small compared to the typical 
wavelength A, one has g [0] ^ — iJA (it is possible, but not necessary at this 
point, to give a precise meaning to the symbol ~). The above linearized equation 
then reduces to a wave equation of speed c, with 

2 , H+H- 

We then use this velocity to nondimensionalize t and ip, 

t = ——, ■>p(X) = — ■ — with Tpo = • 

A/c tpa c 

(the nondimensionalization of ip is obtained via the equation dtip + g' Q = 0). 
We now remark that 



•^''''■-^Qt'^^C,H^]^ and thus ^[CIV'^^c.mi 



'■[C\i^^^Qt\<^H^]^ and thus g [Q^, ^ I^ g ^[eQ^P 
where g"^ and g^ are as defined in (|2.15p and (jl.lSp . Since moreover 



it is straightforward to deduce ()1.19p from ()1.16p . 
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B Proof of ( 12391) . ( 127231) . ( 12:331) . ( I234l) . ( 127351) 
and ( 12:361) 

Before proving these identities, we give some remarks of general interest: 

- The following classical product estimate holds 

Vs > and to > d/2, \fg\„s < |/|^w*o l^l^.. (B.l) 

- For functions defined on the strip S^, the above formula can be used in the 
horizontal directions to obtain 

Vs>0 and to > d/2, \\A'{f9)h < WfhrH^-^^oWgh. (B.2) 

- For all fc e N, we denote by iJ'*''^(5) the space 

k 
3=0 

(B.3) 

- The space H'^^^/^'^{S'^) is continuously embedding in L^H^{S^) (this is a 
variant of the trace lemma) , 

Vs>0, ll/llLs^if=<ll/llH.+i/2,i. (B.4) 

Finally, let us state the following lemma, which will be used several times. 

Lemma B.l. Let to > d/2 and u e H^{S^), g G H''^^{S^y+^ be such that 

\7^'^ ■ P(E±)V^''u = V^* • g, 
"U=o=0' ^nWU=:Fi ="^^■gU==Fl• 

T/len, for alio <s <to + l, 

||A^V^^7.||2 < M||A^g||2 and ||V^*7.|U.a < M{\\A'gh + \\A'-^d,g-e,h). 

Proof. The first estimate is a consequence of Proposition 2.4 of [3 . For the 
second estimate, and since /i^ is assumed to be bounded, it is enough to prove 
that ||A''^^(9^u||2 < M||g||//s.i. This follows immediately from isolating the 
component d^u in the equation V • PV^ u ~ V • g and using the first 
estimate to control the terms involving u. D 

B.l Proof of (\21M . ( 12:231) and ( 12:241) 

B.l.l Proof of dHHD 

The case s = being proved in Prop. 3.4 and Eq. (3.3) of !3|, we focus on the 
case s > 0. Remarking first that |*P'(/j='=|/fs = |*P(A'*-!/'''')|2, we can deduce from 
(Hm with s = that 

Vnm^lH^^ < M||V^*(A^V^)''|| 

< Af(||V^*((A^V*)''~A^0±)|i2 + ||A^V''*(/)±||2, 
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where (A*^='=)'' stands for the solution to (|2.17p with Dirichlet data K''il;^. Re- 
marking that -u = {K''%b^)'^-K''d^ is as in Lcmma lB.ll with g = [A*,P='=]V'' 0=*=, 
we have 



|V'^*((A^V'±)^-A^0±)||2 < ||g||2 



,± 



< M\\^'-^V^' </.±||2, (B.5) 

where we used the commutator estimate (j2.29p to derive the second inequality. 
The end of the proof is then straightforward. 

B.1.2 Proof of dUM]) 

Integrating by parts in (|2.17p one gets easily 

/ A^^^^iA^^a^ / A''(P(E±)V'^*(/)j^)-A''V^*V|, (B.6) 

where 4>i denotes the solution to (|2.17p with Dirichlet data -01, while 

with X 3' smooth compactly supported function equal to 1 in a neighborhood of 
the origin. One then gets 

|(A^^±^i,A^V2)| < |1A^(P(E±)V''*<^±)|12|1A^V''*0||12. 
From (|B.2P one gets, for all < s < io + 1, 

||A^(P(S±)V^*0±)||2 < (l + ||P(S±)-IdLj.^.o+OI|A^V'^*0?||2 
< M||A^V^*0±||2, 

where the second identity follows from (|2.7p : (|2.23p follows therefore from (|2.18p 
and the observation that ||A*V^ "02112 ^ \/Ml*P^2|i/= (Proposition 2.2 of ^3]). 

B.1.3 Proof of dHH]) 

Let us first prove (|2.24l) in the case 1 < s < ip + 1- With the same notations as 
in TO. 1.21 one gets 

([A^^±]0l,A>2) = / [A^P(E±)]V^*Vl•A^V^*0J 

Js± 

+ f P(E±)V''''[A'*0j - (A'^Vi)"] • A'^V^^^J. 
Js± 

Using Cauchy-Schwarz inequality, p.29p (for the first component of the r.h.s.) 
and (IB.SP (for the second component) we get 

|([A^ ^±]0i, A^V2)| < A/IIA^-iV^^^JlUllA^V'^vllls. 
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Since s - 1 > 0, ([^1^ and the observation that jjA^V^ i/;||l2 < y/Jl\^tp2\H- 
yield the resuh. 

Let us now prove (|2.24p for < s < 1. Let us write first that 

([A^ ^ilV-i, A^^2) = -(A^[A, ^±]A-Vi, A^V2) + ([A^+\ ^±]A-Vi, A^V2). 

Since s + f > f , we can use the computations above to show that the second 
term of the r.h.s. satisfies the desired estimate. We thus focus on the first term. 
Proceeding as above, we write 

iA'[A,g^]A-^i:i,A'i;2)= I A"[A,P(E±)]V^*(A-Vi)''-A'^V''*V| 

+ f A'{p{j:^)v''^ [A(A-iVi)'' - V-?]) • A^yf^^i^l 

Since ||A''[A,P(I]±)]|lH-n_>L2(5±) + |lA^(F(I]±)-)||ff=.o^i2(5±) <M (recall that 
< s < 1), we deduce (proceeding as for (jB.SP for the second component) that 



|(A^[A,^±]A-Vi,A^i^2)|<M||A^V^*(A-Vi)''l|2||A^V'^*V4ll2, 
and the result follows as for the case s > 1. 

B.2 Proof of ( 127331) . (12734]) . ( 12735]) and ( 12736]) 

Throughout this proof, we assume that S* is a regularizing diffeomorphism of 
the form S='=(X, z) = {X, z + a^{X, z)) with g^ as in Example WM We denote 
by d?P^{VL) {j e N) the j-th derivative of the mapping ( H> P(S^) in the 
direction h — {hi, . . . ,hj), and by d^ip{h.) the j-th derivative in the direction h 
of C H> 0, where solves (I2.17p . 

B.2.1 Proof of dUMD 

With the same notations as in TO. 1.21 one has 



{A''-^/^g^i;,ip)= f A'*(P(I]±)V''^<?!)±)-A-i/2v 
Js+ 



'%'; 



since ||A ^/^V^ "15^112 ^ ^J'^^^W\2^ it follows after differentiating this identity 
with respect to C, and by a duality argument that for all < s < to + 1 1 one has 

|d^"^±(h)V'|^.-V2 < m'/^^ ||A^d^'iP±(hz)V^*d-'"^0(hzz)||2, 

where the summation is taken over all integers ji and J2 such that ji + J2 = J 
and on all the ji and J2 uplets h/ and h// whose coordinates form a permutation 
of the coordinates of h. The result follows therefore from the estimate, for all 
< s < io + 1/2 (and using Notation I^Jl)) . 

||A«d-'"iF±(h/)V^*d^'^</)(h,,)||2<Me^V7^|/inH=+i/2|h/|H-*o+3/2|Wlff-*o+i/2; 

the rest of this section is thus devoted to the proof of (jB.7p . We first state three 
lemmas and use them to prove the result, and then we prove the lemmas. 

76 



Lemma B.2. For all < s < to + I, k e W , h = {hi, . . . , hk) € i7«v*o+i (!»'*)'= 
and I < I < k, one has 

||A^d'=P±(h)||2 < Me''\hi\Hs+u2\hi\H'-^*o^^- 

Lemma B.3. For all < s < to + 1/2, k e N* , h = {hi, . . . , hk) € iJ*°+2(R'*)'= 
and ip e iJ"'^*o+i(M'^). Then d'=(/)(h) satisfies 

Lemma B.4. Under the assumptions of Lemma \B.3l one also has 

Let us assume first that hi E hj. Then, from (|B.1[) and (IB. 41) . one deduces 



and (jB.7p easily follows from the first two lemmas. 
Now, if hi G h// then we rather write 

||A^d^'ip±(h,)V^*d^'=0(h„)||2 < ||d^'ip±(h)|Uooff.v*o |lA«V''*d^'^0(h„)|l2, 

and Lemma fB. 41 gives (|B.7[) (the bound on ||(i'^P*(h)||ioo^svto being a straight- 
forward consequence of (|B.7p below). 

Proof of Lemma \B.S\ Computing explicitly the matrix P(E ) given by (|2.5p . 
one obtains 

/ Iddxd + d,a^ -V7^Va± \ 
P(S±)= r-,(^^±.T l + M|V^^r ■ (B.7) 

Since C i— > cr* is linear, it is easy to deduce the result from (IB.ip and the 



regularizing properties of a as given in Example 12.41 D 

Proof of Lemma \B.S\ By differentiating ()2.17p k times with respect to C in the 
direction h, one gets 



dV±(h)u.o = 0, a„d'=0±(h)|^^^, = -e, . g|^^^, 



where g is a sum of terms of the form 

d'^-ip±(h/)V^^d'=^</)±(h,,) =: A{\ii,\in), (B.9) 

where fci, ^2 G N, ^2 < k, ki + k2 ~ k, and h/ and h// are respectively fci and 
fc2-uplets whose coordinates form a permutation of the coordinates of h. Since 
by Lemma IB. II we have 

||V^*dV*(h)||H-'0 + V2,l < M|lg|l^w.o + l/2,l, 

77 



we can deduce the result from the fact that 

||A(h/,h//)||^.V*o + l/2,l <M£'=y^|h|jjwto+3/2|Wlff'0 + l, (B.IO) 

which we now turn to prove. From (JB.9P and (jB.2p . one gets 
and therefore, using (jB.71) and the fact that ( i~^ a^ [(] a hnear, 



Since k2 < k, one can deduce (|B.10[) by a simple induction on k (the case k — 
being a consequence of P-lSp ). D 

Proof of Lemma \B.4\ Proceeding as in the proof of Lemma IB. 3| we can check 
that the result follows from 

||A''A(h7,h//)||2 < Me''^\hi\Hs+i/2\hi\Hsv,„+3/2\^iP\H^vt„+i/2, (B.ll) 

that we now turn to prove. We have to distinguish two cases: 
- If hi belongs to h/ we use (|B.2[) and (jB.41) to write 



and (JB.III) follows directly from Lemmas IB. 21 and IB. 31 
- If hi belongs to h//, we rather write 

||A^A(h,,h,,)||2 < ||d^^P±(h,)||i5oH-*o||A^V^*d'=^0(h„)||2, 

and (jB.llI) follows easily by induction as in the proof of Lemma IB. 31 D 

B.2.2 Proof of dHHD 

With the same notations as in i|B.1.2l one has 



(A"-i/2^±V,V5)= / A"+i/2(P(I]±)V^0±)-A-iV 
Js+ 



-iwM^t. 



since \\A ^V^(f^\2 ^ y/JI\'^\2, it follows after differentiating this identity with 
respect to ( and by a duality argument that for all < s < to + 1, one has 

|d^'^±(h)^|^._i/2 < V^^||A^+i/V^P±(h,)V''*d^"^0(h,,)||2, 

where the summation is as in i jB.2.11 The result is therefore a direct consequence 
of (|B.7I) (with a shift of 1/2 derivative). 
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B.2.3 Proof of dUMD 

We proceed as in i )B.1.2l after differentiating (|B.6I) with respect to C- One thus 
gets 

where the summation is the same as in the previous section. Ehiptic estimates 
on (IB.8|) and an induction on j show that 



and the result foUows. 

B.2.4 Proof of pTM]) 

We proceed as for (jOSj) but use ((BJ)) to control \\A'dJ^P^{hj)V>'^ d^^(j)^{hii)\\2. 

C Proof of ( 15341 ) and ( l535l ) 

C.l Proof of (I5A4|) 

We prove ()5.14p in the case i = (that is, J' = 1); the generalization to the 
case i > only requires simple technical modifications and is thus omitted. 
By definition of Ins [U] , one has 

The fact that the second and third term of this identity are controlled by m^(C/) 
is a direct consequence of the last two points of Proposition 14.11 We are thus 
led to control the first term. Let us remark first that, by definition of T [U], 

(E^l • £([Z^KCm),C(o)) = ilV^l • EdZ^lV • (C(a)Z+)),C(a)) 

+ P^H-{lv^j-^ilv^jg{g-r'y-iC(o.)lv^WXia.)) 
= A1 + A2. 

Owing to (|5.6p . the inequality (|5.14p follows from the estimate 

e^fip+p-{\Ai\ + \A2\)<m''{U) 

I + eVp- (1(1 + ^'/'\D\'^')F\l + 1(1 + ^'/^\D\'/^)G\l)] , (C.l) 

with F = C^a)lyj^l and G — C{a)djVlVj^i, that we prove now. 

- Control of A\. We decompose further Ax into A\ = An + A12 + A13 + A14 
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with (summing over the repeated index j, denoting F = (1 + /i^/'*|Z3|^/^)F and 
with g as in (jXTCl) ). 

Ai2 = ti'^^d,[\D\'/\V+]F,Vii'g-^V ■ F), 

A (V+F (^ + V^'-P') ^^^ 0.f r W F^ 



^- - -fe^^' (iv;:^i;.r ) ^[^^^" (1 + .v^i^iv.) - Qp(-)]v ■ ^). 



where *P' = i/imi ^^^d the symbol r is given by (with 5** as in p.3p ) 



^ '^' a + V^^\^\){p+H+s--p-H-s+){l + ^,^/^\^\y^y 

Using Cauchy-Schwarz inequahty, standard commutator estimates and Remark 

12.101 we get that An and A12 are controlled as in (|C.1I) . 

For Ai3, we remark that 0{U) — Ht nii |i/4m i/2-| Tof ^P(''")^"^ ^^ ^ ^^^* order 

operator with skew symmetric principal symbol and that, as a consequence of the 

pseudodifferential estimates of [37], one has /^||0([/) + 0([/)*||i2_>^2 < m^{U). 

It follows that Ai3 is also bounded by the r.h.s. of (|C.1|) . 

For Ai4, one deduces from Cauchy-Schwarz inequality that 

\A^,\ < \VlFU{\[^'d-' (^ ^ ^m\D\i/2) ~ OpiT)]V ■ F\^ 

and Corollarv l3.12l (with /c = for the first term and k ~ 1 for the second one) 
implies that {A^l < -m^{U)\F\2, from which (jC.ip follows for for A14 and 
therefore for Ai. 
- Control of A^- Let us first rewrite A^ as (with F = C(q)[11*]) 

Using ([Z;25)) with u = (^"^V • F and v = |F=^], and using Remark [2?T0l vields. 
with F as above, 

|^2| < ^m^(^)|F|i 
which implies (|C.ip . 
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C.2 Proof of dmi]) 

As for (|5.14p . we just consider the case i = (i.e. J'' — Id). We first write 

M 
= B1 + B2. 

For Bi, we use the fact that g = -^Q^ ° •^[C] "^ (see Remark [2. 22p to get 

Si = --(Z+-VJ[C]g,^+.9), 
with5 = J[C]"V(a)- Since J[C] =/?+ Id - p"^(^")"i^ + , we deduce that 

with 5 = {g^)^^g^g. It follows therefore from (|2.25p . Lemma [2. 191 and Propo- 
sition H^ that |Si| < m^{U). For B2, we just write 



-(K±l-Vg«,C^-g«) 



^2 

M 

with 5' = (^-)-i^?/;(„). Using ((2:25)) again, we also get IB2I < m^{U), and 
the proof of (|5.15l) is complete. 
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